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Sobald in den Koeffizienten g(x) und o(x) eines Eigenwertproblems 
d df 
= Hs Ein] 0 
Zee CAR 26 G)) f | (1) 


ein entsprechender Parameter auftritt, so kann es, wie in §2 und § 3 angegeben wird, 
in ein anderes Eigenwertproblem „umgeordnet“ werden, das die gleichen Eigenfunktionen 
besitzt, wie das ursprüngliche (1), falls es nur in dem gleichen Grundgebiet und für die 
gleichen Randbedingungen formuliert wird. Die mit Hilfe der Polynommethode lösbaren 
Eigenwertprobleme ergeben nicht nur lineare umgeordnete Eigenwertprobleme von der | 
Gestalt (1), sondern auch quadratische von der Gestalt 


d df 4 ar 20 
Can ee A? o” (x)) f = 0. 


Aus den linearen umgeordneten Eigenwertproblemen kann man für die Eigenlösungen 
des urspriinglichen Eigenwertproblems neue Orthogonalitatsrelationén von der gewohnten 


Gestalt 


Xs : 
fra | ed = o (=m, 5 
21 = a 
edoch mit anderen Gewichtsfunktionen erhalten. 
Wenn ein gegebenes Eigenwertproblem in ein lineares umgeordnet werden kann, 
können seine Eigenfunktionen auch als Lösungen eines linearen zweiparametrigen Eigen- 


wertproblems 
d d 
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angesehen werden ($ 4). Bei den quadratischen und zweiparametrigen Eigenwertproblemen 
treten in den Orthogonalitätsrelationen die Eigenwerte der entsprechenden Eigenfunk- 
tionen auf. Dies hat aber zur Folge, dass die Koeffizienten der eventuell vorhandenen 
Entwicklungen willkürlicher Funktionen nach Eigenfunktionen aus unendlich vielen 
linearen Gleichungen mit unendlich vielen Unbekannten zu berechnen sind. Aus den 
Eigenlösungen zweiparametriger Eigenwertprobleme kann man jedoch in mannigfacher 
Weise Sätze von unendlich oder endlich vielen Eigenfunktionen auswählen, die neue 
Orthogonalitätsrelationen von der gewohnten Gestalt (2) liefern. 

In $5 wird gezeigt, dass vom mathematischen Standpunkt der in der Faktorisie- 
rungsmethode auftretende Parameter m als die „Eigenwert-Quantenzahl“ eines umgeord- 
neten linearen oder quadratischen Eigenwertproblems anzusehen ist. Ferner wird 
darauf hingewiesen, in welcher Richtung man eine Erklärung der merkwürdigen Tat- 
sache zu suchen hat, dass in den faktorisierbaren Eigenwertproblemen die Eigenwerte 
für um ganze Zahlen sich ändernde Parameter m stets der gleichen Folge von Zahlen 
zu entnehmen sind. 


$ 1. Eigenwertprobleme der Quantentheorie, deren Lösungen mit Hilfe der Polynom- 
methode herstellbar sind 


Im folgenden beschäftigen wir uns mit Eigenwertproblemen, der Quantentheorie, 
die durch Differentialgleichungen zweiter Ordnung definiert und mit Hilfe der Sommer- 
feldschen Polynommethode (Sommerfeld 1939, vgl. auch Rubinowicz 1949, 1950, 
1959,M) lösbar sind. Wir können sie stets in der selbstadjungierten Normalform 


^ وام‎ FE — (g(a) — den) £ = 0 a) 


voraussetzen. Wir nehmen an, dass der Grundbereich des Eizenwertproblems durch 
zwei singuläre Punkte x, und x, der Differentialgleichung (1.1) begrenzt wird, wie dies 
in der Quantentheorie in der Regel der Fall ist. Wir beschäftigen uns im folgenden 
ausschliesslich mit dem Fall diskreter Eigenwertspektren. Bezüglich des Verhaltens 
der Eigenlósungen in den Begrenzungspunkten x, und x, des Grundbereiches setzen 
wir daher im folgenden voraus, dass hier ein eventuelles Unendlichwerden der Eigen- 
lösungen die Konvergenz der Integrale der Normierungs- und Orthogonalitätsbedingung 


f f 0 ds = s, (13) 


nicht behindern darf. Da zur Ableitung von (1.2) die Forderung notwendig ist, dass 


lim G, = lim 6, f (1.2a) 


X>% XX: 


ist, wo Gin durch 


Cnn = Pp (u — fa 2 
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definiert wird, ist somit die Bedingung (1.2) den Grenzbedingungen (1.2a), (1. 2b) i in 
den Endpunkten x, und x, des Grundbereiches äquivalent. Für die weiteren Über- 
legungen mag ausdrücklich darauf hingewiesen werden, dass in den Grenzbedingungen 
(1,2a), (1.2b) nur der Koeffizient p(x) des Eigenwertproblems (1.1) auftritt. Sie hängen 
daher weder von dem Eigenwertparameter A. noch von irgendwelchen in den Koeffi- 
zienten q(x) und g(x) auftretenden Parametern ab. Wie sich in § 2 und § 3 bzw. in $4 
herausstellen wird, können daher die Grenzbedingungen (1.2a), (1.2b) formell auch für 
die (linearen und quadratischen) umgeordneten bzw. für die mehrparametrigen Eigen- 
wertprobleme unverändert beibehalten werden. Hingegen werden die Normierungs- 
und Orthogonalitätsbedingungen (1.2) in den verschiedenen angeführten Fällen 
verschiedene Gestalten annehmen. 

Um für die Funktion q(x) eine kurze Bezeichnung zur Verfügung zu haben, 
wollen wir sie in der Folge als die Potentialfunktion bezeichnen, da sie in sehr vielen 
Fällen mit dem im Hamiltonschen Operator auftretenden Potential im Zusammen- 
hang steht. 

Als Beispjele solcher Eigenwertprobleme führen wir an: das Eigenwertproblem 
der zugeordneten Kugelfunktionen 


d 1 df m? i 
snene [1m \r=o (1.3) 


sowie das der Radialfunktionen des Einelektronenatoms 


oe ا‎ 0. (1.4) 


Der Grundbereich des ersten Eigenwertproblems ist (—1,+1) und der des zweiten 
(0,20). Wie der Vergleich von (1.3) bzw. (1.4) mit (1.1) lehrt, haben die Koeffizienten 
,ووم‎ o im Falle des Eigenwertproblems der zugeordneten Kugelfunktionen die Gestalt 


m ? 


e =1 (1.5)‏ ,= و «ه-1 - )م 


und werden im Falle des Eigenwertproblems der Radialfunktionen durch 


p) = x°, واو‎ - 1) +1-2 g(x) = x? (1.6) 


gegeben. 
In der Quantentheorie wird vorausgesetzt, dass die in der N q(x) 
in den beiden Fällen auftretenden Parameter m und l ganze Zahlen sind. Um jedoch . 
diese Eigenwertprobleme mit Hilfe der Polynommethode aufzulösen, ist diese Vor- 
aussetzung nicht notwendig und würde uns auch bei den folgenden ‚Überlegungen 
behindern. Wir wollen daher nur annehmen, dass m reell und / positiv ist, so dass 
die Parameter m? und i(1--1), die in den Potentialfunktionen g(x) in (1.5) und (1.6). 
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auftreten, positive Werte haben. Die Eigenfunktionen des Eigenwertproblems (1.3) 
werden für positive m-Werte durch die zugeordneten Kugelfunktionen 


m (21)! d m/2 p loq ay 2321 top RS 
P; (x) = SH (L — m)! x! (1 = x?) /2 of = 9 (i m), 9 (i m 1); l + 
+ E : ipi) (1.7) 
und die Eigenwerte durch 
۸= ¿(i+ 1) (1.8) 


gegeben. Damit die Kugelfunktion (1.7) eine Lösung des Eigenwertproblems der 
zugeordneten Kugelfunktionen (1.3) ist, muss die in (1.7) auftretende hypergeometri- 
sche Funktion ein Polynom sein. Es muss daher angenommen werden, dass 


O O OG SO Ec (1.9) 


ist. Die Differenz der beiden Zahlen 7 und m muss somit gleich sein einer nicht 
negativen ganzen Zahl n. Da n darüber entscheidet, von welchem Grade das in 
der Eigenfunktion (1.7) auftretende und aus der hypergeometrischen Funktion 
oF (a,b; c;£) entstehende Polynom ist, wollen wir n als die Polynomquantenzahl 
bezeichnen. Es sei hier betont, dass nur Kugelfunktionen deren Parameter / und m 
der Bedingung (1.9) genügen, deren Parameterdifferenz l—m also durch ein eganze 
Zahl n gegeben ume als Lósungen des Eigenwertproblems (1.3) angesehen wa 
können. 

Im Falle des Eigenwertproblems (1.4) der Radialfunktionen werden die Eigen- 
funktionen des diskreten Eigenwertspektrums bestimmt durch 


R. (z) = xt el F (— (n — I— 1);21+2;%x), x = 2/na,, (1.10) 
wo 
n = n, + l4 1. (1.11) 


n, = n—l—1 hat hier die Rolle der Polynomquantenzahl inne. Die Differenz der 
beiden Parameter n und J muss also eine ganze Zahl sein. Die Eigenwerte werden durch ` 
ipee (1.12 
i223 DE e 

gegeben. 
Als ein charakteristischer Zug der Eigenwertprobleme, die mit Hilfe analytischer 
Methoden, insbesondere also der Polynommethode lösbar sind, ist die in unseren 
beiden Beispielen festgestellte Tatsache anzusehen, dass ihre Lösbarkeit erhalten — 


bleibt, wenn wir die in der Potentialfunktion g(x) oder auch in den Koeffizienten p(x) 


und o(x) auftretenden Konstanten in gewissen Bereichen beliebig ündern. Die einzige 
in den Lösungen auftretende Konstante, die unbedingt ganzzahlig sein muss, ist 
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die Polynomquantenzahl, die die gewöhnlichen oder konfluenten hypergeometrischen 
Funktionen „F(a, b; c;£) bzw. ,F,(a; c;£) im Falle eines diskreten Eigenwertspektrums 
in ein Polynom verwandelt, wenn a oder b bzw. a gleich einer negativen ganzen Zahl 
gesetzt wird. Die Quantenzahlen, die die Eigenwerte festlegen (in unseren beiden 
Beispielen / bzw. n) und die wir daher als die Eigenwertquantenzahlen bezeichnen 
wollen, müssen hingegen nicht ganzzahlig sein damit zur Lösung der betreffenden 
Eigenwertprobleme die Polynommethode anwendbar ist. In der Quantentheorie 
werden jedoch die Eigenwertquantenzahlen in der Regel durch ganze Zahlen gegeben. 


82. Umordnung von Eigenwertproblemen 


Aus den mittels der Polynommethode lósbaren Eigenwertproblemen (1.1) 
kann man im allgemeinen andere, ebenfalls mit Hilfe dieser Methode lósbare Eigen- 
wertprobleme herleiten. Es muss nur in der Potentialfunktion q(x) und eventuell 
auch in der Gewichtsfunktion o(x) ein in p(x) nicht enthaltener Parameter, er heisse 
m, derart auftreten dass die Beziehung 


q(x) — Ae(x) = q'(x) — A'o'(x) 


gilt. Dabei ist A’ = A'(m) eine Funktion des Parameters m, von dem jedoch g’(x) 
und o'(x) nicht abhängen. Wir können dann das gegebene Eigenwertproblem (1.1) in 
das Eigenwertproblem!) 


= p) E — (aa) — Nel f = 0. 


q'(x) = q(x) — Ae(x) + A'(m) e'(x) (2.1) 


„umordnen“, wobei der Parameter m weder in der neuen Potentialfunktion qg'(x) 
noch in der neuen Gewichtsfunktion o'(x) enthalten ist. 


Nehmen wir nun an, dass das ursprünglich gegebene Eigenwertproblem (1.1) für 
einen gewissen Wert des Parameters m = m, gelöst wurde und man dabei die zum 
Eigenwert A = A, gehörige Eigenfunktion f = f, erhielt. Setzen wir voraus, dass in 
den beiden Eigenwertproblemen, dem ursprünglichen und dem uingeordneten, das 
Grundgebiet und die Randbedingungen (1.23), (1.2b) die gleichen sind. Nehmen 
wir ferner an, dass wir das umgeordnete Eigenwertproblem (2.1) für den Fall lösen, 
wo A = A, ist, so ist A’ = A'(m,) ein Eigenwert des umgeordneten Eigenwertproblems 
(2.1) und die zu diesem Eigenwert gehörige Eigenlósung des umgeordneten Eigen- 
wertproblems (2.1) wird durch die zum Eigenwert A — A, gehörige ê 


. f = f, des ursprünglichen Eigenwertproblems (1.1) gegeben. 


1) Ich habe mir erlaubt die Bezeichnung „umgeordnete Eigenwertprobleme“ einzuführen, weil 


. ich auf diese Problemgruppe in der mir zugänglichen mathematischen Literatur keine Hinweise und 
— daher auch keinen Namen gefunden habe. 
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Beachtet man, dass mit Hilfe der Polynommethode das ursprüngliche Eigen- 
wertproblem (1.1) für beliebige, in einem gewissen Bereiche gelegene Werte des Para- 
meters m und daher auch von A’(m) und das umgeordnete Eigenwertproblem (2.1) für 
beliebige, in einem gewissen Bereiche gelegene Werte von A lösbar ist, so kann man 
sagen: Die Gesamtheit der Eigenfunktionen des ursprünglichen Eigenwertproblems 
(1.1) (für beliebige Werte des Parameters 4’) und des umgeordneten (2.1) (für beliebige 
Werte des Parameters A) ist die gleiche, gleiches Grundgebiet und gleiche Rand- 
bedingungen (1.2a), (1.2b) vorausgesetzt. 

Man kann ferner zu jedem gegebenen ursprünglichen Eigenwertproblem ein 
bestimmtes umgeordnetes angeben, das die gleiche Eigenfunktion, wie das ursprüng- 
liche besitzt. Umgekehrt kann man zu jedem gegebenen umgeordneten Eigenwert- 
problem ein entsprechendes ursprüngliches mit der gleichen Eigenfunktion finden. 

Die Eigenwertprobleme (1.1) und (2.1) werden wir in der Folge als lineare Eigen- 
wertprobleme bezeichnen, da in ihnen die Eigenwertparameter nur in der ersten Potenz 
auftreten. 

Im allgemeinen kann man nicht erwarten, dass in dem umgeordneten Eigen- 
wertproblem (2.1) die Gewichtsfunktion o’(x) alle Bedingungen erfüllt, die wir ihr 
sonst auferlegen. Sie kann im Grundgebiet ihr Vorzeichen ändern, so dass z.B. nicht 
alle Normierungsintegrale ein positives Vorzeichen haben. Oder sie kann eine Singula- 
rität aufweisen, die zur Felge hat, dass alle oder auch nur einige Integrale in der 
Normierungs- und Orthogonalitätsbedingung für die Eigenfunktionen nicht konver- 
gieren. Oft begegnet man auch dem Falle, dass das System der Eigenfunktionen 
nicht vollständig ist, weil z.B. nur eine endliche Anzahl von Eigenfunktionen auftritt. 
Jeder dieser Umstände bedingt jedoch, dass den entsprechenden umgeordneten Eigen- 
wertproblemen keine physikalische Bedeutung zugeschrieben werden kann. 

Die beiden Eigenwertprobleme (1.3) und (1.4) können in der nachstehenden 
Weise in lineare Eigenwertprobleme umgeordnet werden. Substituiert man in (1.3) 
für A den Wert (1.8) und setzt ferner 


A censet (2.2) 


(Minuszeichen damit ọ'(x)>0 und daher das Normierungsintegral positiv ist!), so 
erhält man aus dem Eigenwertproblem (1.3) der zugeordneten Kugelfunktionen das 
umgeordnete Eigenwertproblem 


d df 1 
——— 2 m Fe وملام شع‎ LÀ s=. 
i ler: | 11+1) (=٥ (2.3) 
Ebenso ergibt sich aus (1.4) mit Rücksicht auf (1.12), wenn man 
m EA لو ف‎ )2.4( 
setzt, das umgeordnete Eigenwertproblem x 


d df, x x | 
mdi peus SF C E a = 
dx” dx | = a?n? x) pao 25 


< 
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Nehmen wir an, dass die beiden umgeordneten Eigenwertprobleme (2.3) und 
(2.5) in dem Grundgebiet (—1, +1) bzw. (0,00) der ursprünglichen Eigenwert- 
probleme (1.3) bzw. (1.4) definiert sind und auch den gleichen Randbedingungen 
(1.2a), (1.2b) genügen, so werden die Eigenfunktionen der beiden umgeordneten 
Eigenwertprobleme durch die Eigenfunktionen (1.7) bzw. (1.10) der ursprünglichen 
Eigenwertprobleme gegeben. Die Eigenwerte werden durch (2:2) bzw. (2.4) bestimmt. 
wobei m in (2.2) gemäss (1.9) durch m = I— n und / in (2.4) gemäss (1.10) durch 
1l n—n,— 1 gegeben ist. 

Umgeordnete Eigenwertprobleme sind als solche für den Physiker nur dann 
von Interesse, wenn sie eine physikalische Bedeutung haben. Dies trifft jedoch meist 
nicht zu, so dass man vielleicht geneigt ist die obigen Überlegungen als eine müssige 
Spielerei anzusehen. Dies ist jedoch keineswegs der Fall, insbesondere wenn das 
umgeordnete Eigenwertproblem linear ist. Man kann dann nümlich aus dem umge- 
ordneten Eigenwertproblem neue Orthogonalitätsrelationen von der üblichen Gestalt 
(1.2) für die Eigenfunktionen des ursprünglichen Eigenwertproblems herleiten, die 
sich eventuell als nützlich erweisen können. Und zwar entspricht jeder Möglichkeit 
aus einem gegebenen Eigenwertproblem ein umgeordnetes lineares herzustellen eine 
neue Orthogonalitätsrelation von der Gestalt (1.2). 

Als erstes Beispiel führen wir die Orthogonalitäts- und Normierungsbedingungen 
für das ursprüngliche und das umgeordnete Eigenwertproblem der zugeordneten 
Kugelfunktionen an 
Ti MET 
| Pray Pr 6) dx = by (1— V = ganze Zahl), (2.62) 
-1 


+1 


f Pra) PI (2) + 


-1 


i = dx = Onn’ (m — m’ = ganze Zahl). (2.6b) 


Da gemäss (1.9) die Differenz Í — m eine ganze Zahl sein muss, treten, wie leicht er- 
sichtlich, in der Orthogonalitütsbedingung (2.6a) bzw. (2.6b) nur zugeordnete Kugel- 
funktionen auf, deren untere bzw. obere Indizes J, P bzw. m, m’ sich voneinander 
nur um ganze Zahlen unterscheiden. Man bemerkt auch, dass die Normierung der 
zugeordneten Kugelfunktionen in den beiden Fällen (2.62) und (2.6b) eine verschie- 
dene ist. Jede der beiden Orthogonalitatsrelationen (2.62) und (2.6b) gilt jedoch selbst- 
verständlich für beliebig normierte Kugelfunktionen. 

Die Orthogonalitäts- und Normierungsbedingungen für das ursprüngliche und 
das umgeordnete Eigenwertproblem der Radialfunktionen des Einelektronenatoms 


sind durch 


8 A 
f Ra (x) R ,, (x) x? dx = Ony  (n— n' = ganze Zahl), (2.7a) 
E 


f Ru) Ry وما‎ d< = ó, (0— I” = ganze Zahl) - (2.7b) 


- 
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festgelegt. Mit Rücksicht darauf, dassnach (1.11) die Differenz n — l eine ganze Zahl 
sein muss, kann man leicht schliessen, dass die Differenzen n — n' in (2.7a) und 1 — I’ 
in (2.7b) ebenfalls ganze Zahlen sein müssen. 

Die betrachteten Eigenwertprobleme sind Beispiele für den Fall, wo die Anzahl 
der Eigenfunktionen im umgeordneten Eigenwertproblem eine endliche ist. Bei 
gegebenem A, d.h. im Falle der zugeordneten Kugelfunktionen bzw. im Falle der 
Radialfunktionen bei festgelegtem / bzw. n erhalten wir nämlich eine endliche An- 
zahl von Eigenfunktionen, da ja |m| >> I bzw. 1 <n sein muss. 

Ein solcher Tatbestand scheint in sehr vielen Fällen aufzutreten. Eigenwert- 
probleme, die durch Umordnung aus Eigenwertproblemen mit einem vollständigen 
System von Eigenfunktionen entstehen, besitzen im allgemeinen nur eine endliche 
Anzabl von Eigenfunktionen. Dies wird durch den in $ 5 besprochenen Zusammen- 
hang zwischen der Polynomquantenzahl und den Eigenwertquantenzahlen des ur- 
sprünglichen und des umgeordneten Eigenwertproblems bedingt. 

Es gibt auch Eigenwertprobleme, die mit Hilfe der Polynommethode lösbar 
sind und die in mehrere verschiedene Eigenwertprobleme umgeordnet werden können. 
Ein Beispiel dafür sind die verallgemeinerten zugeordneten Kugelfunktionen (vgl. 
Bateman 1949, S. 389), die neuerdings von Kuipers und Meulenbeld (1957, vgl. auch 
Kuipers 1958, Meulenbeld 1958, 1959) näher untersucht wurden. Sie sind Lösungen 
der Differentialgleichung 


NR E TE stern | حم‎ (2.8) 


die im Falle m? = n? in die Differentialgleichung (1.3) der gewöhnlichen zugeordneten 
Kugelfunktionen übergeht. Je nachdem 

(a) A = ¿(i + 1) oder 

(b) A’ = — 4 m? oder 

() A= — pn? 


gesetzt wird, erhalten wir drei verschiedene Eigenwertprobleme?) von denen jedes 
als das ursprüngliche und die beiden restlichen als die umgeordneten angesehen 
werden können. E 

Eine Lósung der Differentialgleichung (2.8) der verallgemeinerten zugeordneten 


Kugelfuriktionen, die im Bereiche (— 1, +1) reell ist, wird durch (vgl. Bateman 1949, 
S. 389) 


pron (x) = (1 = AA + a) (1 +} (m + n), 14 Y (m + n) + 
+l1;m+1;3(1—-%) (2.9) 


gegeben. Damit vu als eine Eigenfunktion eines der drei. Eigenwertprobleme (a), 
2) Die Eigenwertprobleme (b) und (c) gehen TM fibers falls m* mit n? vertauscht und gleich- 
zeitig das Vorzeichen von x geändert wird. 
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(b) Mr (c) angesehen werden kann, muss die hier auftretende hypergeometrische 
Funktion in ein Polynom übergehen und daher im Falle positiver /-, m- und n-Werte 


l— + (m +n) ع‎ 0,12, ... (2.10) 
sein. Die Parameter /, m, n müssen zwar die Bedingung (2.10) erfüllen, jedoch nicht 
ganzzahlig sein. 

Im Falle n= m gehen die verallgemeinerten zugeordneten Kugelfunktionen 
(2.9) in die gewóhnlichen zugeordneten Kugelfunktionen P7‘(x) über, da die letzteren 
nicht nur durch (1.7) sondern (in Darwinscher Normierung) auch durch 

(m + D! : 

Pr. (x) = 2m m1(l— m)! (1— x?) F; (m — D m — l + il m — 1;3 (1— x)) 

darstellbar sind. 


Die den drei Eigenwertproblemen (a), (b) und (c) entsprechenden Orthogonali- 
tátsrelationen lauten 


+1 
(a) [P= O) Pp" @) dx 0, (ee 
ze 


+1 
(b) | کح منرم مسجم‎ = o, (m — m' = + 2, + 4, + 6, ...) 
= 
+1 
n m,n’ dx f i 
(c) fr (x) Pr» ا‎ e (n—n = +2, £4, +6, ...) 


-1 
Die Bedingungen für die Quantenzahlendifferenzen l, m, n ergeben sich dabei aus 
(2.10). 

Da in den mit Hilfe der Polynommethode herstellbaren Eigenlösungen höch- 
stens drei willkürliche Parameter auftreten, (vgl. Rubinowicz M), so kann man daraus 
schliessen, dass sich aus solchen Eigenwertproblemen höchstens zwei umgeordnete 
Eigenwertprobleme herstellen lassen. 

Die Bezeichnung ursprüngliches und umgeordnetes Eigenwertproblem zeichnet 
keines von ihnen in irgend einer Weise aus. Man kann ja stets ein umgeordnetes 
Eigenwertproblem als das ursprüngliche und ein ursprüngliches als das umgeordnete 
Eigenwertproblem ansehen. 


5 3. Quadratische Eigenwertprobleme 


Nicht in jedem Spezialfalle der mit Hilfe der Polynommethode lösbaren Eigen- 
wertprobleme tritt ein Parameter m in dem Ausdruck q(x) — Ao (x) des ursprüngli- 
chen Eigenwertproblems in der Weise auf, dass es sich in ein lineares Eigenwert- 
problem (2.1) umordnen lässt. Unter den mit Hilfe der Polynommethode lósbaren 
Eigenwertproblemen sind auch solche enthalten, die sich nur in Eigenwertprobleme 
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von der Gestalt 
^ pa) É — 0) — e a) — 270" (9) f = 0, 
q (6) = qla) — te (a) + 3' “و‎ (x) + 229” (a) (3.1) 


umordnen lassen, in denen also der Eigenwertparameter A’ quadratisch auftritt und 
die wir daher kurz als quadratische Eigenwertprobleme bezeichnen wollen. Von m ist 
dabei nur der Eigenwertparameter A’ abhängig, während p, g', o', o" von m unab- 
hängig sind. 

Selbstverständlich muss man voraussetzen, dass o (x) und o" (x) nicht zueinander 
proportional sind: o'(x) = c o" (x). Im Falle o'(x) = c ọ''(x) kann man ja cd’ +4’? 
als einen Eigenwertparameter auffassen und erhält so ein lineares Eigenwertproblem 
mit der Gewichtsfunktion g” (x). 

Zunüchst wollen wir uns überzeugen, dass ein quadratisches Eigenwertproblem 
(3.1) nur unter besonderen Voraussetzungen nur reelle Eigenwerte 4, besitzt, trotz- 
dem es durch eine selbstadjungierte Differentialgleichung gegeben wird. Um diese 
Bedingungen abzuleiten, bezeichnen wir mit f,(x) die zum Eigenwert A, gehörige 
Eigenlösung von (3.1). Multiplizieren wir (3.1) für f= f, und A’ = A, mit f; und 
subtrahieren davon den mit f, multiplizierten konjugiert komplexen Wert von (3.1) 
für f = f, und X = 44, so erhalten wir 


dex 4 Ps, D 4 "T * 
d; Cm = a) + A — An) e" (1, f (3.2) 
wo in Übereinstimmung mit (1.2b) 
7 «df, df, 
Bee 4£ ر‎ 
ist. Unter der (1.2a) entsprechenden Voraussetzung 


lim G,, = lim G,, 


xx, XX 


ergibt dann eine Integration von (3.2) über das Grundgebiet (x,, x) die Beziehung 


MA) f fs f. le (®) ++ 27) 0" (9)] dx = 0. (3.3) 


Um zu beweisen, dass A, reell d. h. A, = A, ist, darf das Integral in (3.3) für den be- 
treffenden n-Wert nicht verschwinden. Dies tritt Sicherlich ein, wenn die eckige 
Klammer in (3.3) im Grundbereiche nicht verschwindet also z. B. wenn Q'(x) und 
o” (z) innerhalb des Grundbereiches gleiche (bzw. verschiedene) Vorzeichen haben 
und A, + A4 für das betreffende n nicht negativ (bzw. nicht positiv) ist. 

Und nun «oll festgestellt werden, in welchem Sinne man von einer Orthogonalität 
zweier zu verschiedenen Eigenwerten gehörigen Eigenfunktionen sprechen kann. 
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Um die Orthogonalitätsbedingung abzuleiten, multiplizieren wir (3.1) für f — 
und 4 = À, mit f, und subtrahieren davon den konjugiert komplexen Wert der 
mit f, multiplizierten Gleichung (3.1) für J f, und 2' = A,. Beschränken wir 
uns auf den Fall reeller Eigenwerte, so erhalten wir 


de om — [Ag A1) e (n) AP OA (3.4) 


wo ebenso wie in (1.2b) 


ist. Integration der Beziehung (3.4) über den Grundbereich ergibt dann mit Rück- 
sicht auf 4, Z 27, die Orthogonalitätsbedingung 


OO. e UL +4) 0” (QE dx 0, (3.5) 


vorausgesetzt, dass die Randbedingung (1.2a) 


lim G,,, = lim Gn 
xx X + Xs 
srfüllt ist und daher das Integral (3.5) konvergiert. Auffallend an der Orthogonalitäts- 
bedingung (3.5) ist ihre Abhängigkeit von den Eigenweriparametern der beiden hier 
ıuftretenden Eigenfunktionen. 
Diese Abhängigkeit bedingt jedoch einen schweren Nachteil der quadratischen 
Eigenwertprobleme gegenüber den linearen. Selbst wenn das quadratische Eigen- 
wertproblem (3.1) ein vollständiges System von Eigenfunktionen besitzt und daher 


jede hinreichend reguläre Funktion f (x) in der Gestalt 
Fl) = Dita Ca 


darstellbar ist, so verhindert die Abhängigkeit der Orthogonalitätsbedingung (3.5) 
von den Eigenwertparametern eine einfache Berechnung der Entwieklungskoeffizienten 
c„. Aus (3.5) kann man nämlich für die einzelnen Koeffizienten c, keine einfachen 
Ausdrücke, wie im Falle der linearen eigenwertprobleme erhalten. Zur Berechnung 
der c, ergibt sich nämlich aus (3.5) ein unendliches System von linearen Gleichungen, 
von denen jede im allgemeinen unendlich viele Koeffizienten c, enthält. 

Als ein Beispiel eines mit Hilfe der Polynommethode lösbaren Eigenwertproblems, 
das in ein quadratisches umgeordnet werden kann, möge das beim symmetrischen 


Kreisel auftretende Eigenwertproblem 


d 6  fe+r dr: |- | 
2-32-| 4 x(1 — x) As 3 0 
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dienen. Setzt man hier für den Eigenwertparameter A seinen Wert 


T wenn T >T; 
, 
T wennt >T 


KEN يي‎ lan 


ein, so kann man in (3.6) entweder v oder 7' als den Eigenwertparameter A’ des umge- ' 
ordneten Eigenwertproblems ansehen, das offenbar quadratisch ist. 

Eigenwertprobleme, in denen der Eigenwertparameter in einem höheren als dem 
zweiten Grade auftritt, können wir ausser Betracht lassen. Es sind nämlich keine 
Beispiele bekannt, wo sie aus den mit Hilfe der Polynommethode lösbaren linearen 
Eigenwertproblemen durch Umordnung herstellbar sind. 

Die Möglichkeit, Eigenwertprobleme, die mittels der Polynommethode lösbar 
sind, umordnen zu können, ist für das Verständnis der Faktorisierungsmethode von 
Bedeutung, wie wir dies in §5 näher begründen werden. 


§ 4. Zweiparametrige Eigenwertprobleme 


Wir wollen uns nun mit einer Gruppe von Eigenwertproblemen befassen, die 
für uns deshalb von Interesse sind, weil sie aus allen Eigenwertproblemen gewonnen 
werden können, die sich umordnen lassen. Sie werden passend als mehrparametrige 
Eigenwertprobleme bezeichnet, da sie sich von den einparametrigen Eigenwertpro- 
blemen (1.1) dadurch unterscheiden, dass in ihnen mehrere Eigenwertparameter 
auftreten. Insbesondere interessieren wir uns für die linearen, zweiparametrigen 
Eigenwertprobleme, die aus den einparametrigen hergestellt werden können, sobald 
sich diese in lineare Eigenwertprobleme umordnen lassen. Wir können sie stets in der 
selbstadjungierten Normalform 


RR (A) 


hinschreiben. Wir wollen annehmen, dass das Eigenwertproblem (4.1) in einem durch 
zwei singuláre Punkte x, und x, dieser Differentialgleichung begrenzten Grundgebiet 
(wi, xə) gegeben ist. In den Endpunkten des Grundgebietes sollen sich die Eigen- 
funktionen so verhalten, wie dies den Bedingungen (1.2a), (1.2b) d. h. den beiden 
Eigenwertproblemen entspricht, die wir aus (4.1) erhalten, wenn wir hier entweder 
AY oder 2 als gegebene Konstante ansehen. Wir stellen uns die Aufgabe, zusam- 
mengehörige Paare der beiden Eigenwertparameter A) und A(? zu finden, für die 
sich Eigenlósungen, d. h. innerhalb des Grundgebietes (x,, x.) reguläre Funktionen 
angeben lassen, für die die Integrale einer der beiden angegebenen Normierungs- 
und Orthogonalititsbedingungen konvergieren. ; 

Wir nehmen dabei an, dass ähnlich wie beim quadratischen Eigenwertproblem 
(3.1) die beiden Gewichtsfunktionen ọ,(x) und g(x) nicht zueinander proportional 
sind, d. h. o,(x) Z c os(x) ist. Im Falle o,(x) = c 0,(x) geht ja das zweiparametrige 
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Eigenwertproblem (4.1) in ein einparametriges (1.1) mit dem Eigenwertparameter 
c AU + A9 und der Gewichtsfunktion Q(x) über. 

Eine Eigenlösung f, des Eigenwertproblems (4.1) gehört stets zu einem Paar 
zusammengehöriger Eigenwerte AU und AO. Zur Bezeichnung der Eigenfunktionen 
und ihrer beiden Eigenwerte verwenden wir im allgemeinen einen und denselben 
Index z. B. n oder m. 

Einparametrige lineare Eigenwertprobleme, die durch selbstadjungierte Opera- 
toren gegeben werden, haben stets reelle Eigenwerte. Bei zweiparametrigen Eigen- 
wertproblemen solcher Operatoren ist das nicht immer der Fall. Um dies einzusehen, 
multiplizieren wir (4.1) für f = f, und AP = 4M, 49 — 4 mit f? und subtrahieren 
davon den mit f, multiplizierten konjugiert komplexen Wert dieses Ausdruckes. 
Wir erhalten dann 


E Cm = [0P — 20) وين‎ + GP =I") eA. CS 


wo im Einklang mit (1.2b) 


Gun D (^ E — fn 7 


ist. Integrieren wir (4.2) über das Grundgebiet (x1, x), so ergibt sich 
OP — AD) f £f, ex aad OS — 12) J fen) de 0, 5 (3) 


sobald die Randbedingung (1.2a) 
lim G,,, = lim G,, 


XX E27 
erfüllt ist, so dass die in (4.3) auftretenden Integrale konvergieren. 

Da die Integrale in (4.3) reell sind, stellt (4.3) eine Beziehung zwischen den Ima- 
ginärteilen 1 (AD — 4°) bzw. 4 (4% — 22") der beiden zusammengehörigen Eigen- 
werte A‘) und 4% dar. Setzen wir voraus, dass 0,(x) und g(x) und daher auch die 
beiden in (4.3) auftretenden Integrale positiv sind, so müssen die Imaginärteile von 
AD und 42) verschiedene Vorzeichen haben. Ferner ergibt (4.3) die sehr wichtige 
Tatsache, dass sobald unter der Voraussetzung 0,(x) > 0, os(x) > 0 einer der beiden 
zusammengehörigen Eigenwerte A(? oder 4% reell ist, es auch der andere sein 
muss. 

Falls wir nicht ausdrücklich das Gegenteil betonen, wollen wir im folgenden 
nur Eigenfunktionen betrachten, die zu reellen Eigenwerten 4% und A gehören. 
Nur solche kónnen ja im allgemeinen für physikalische Überlegungen interes- 
sant sein. | | 

Nehmen wir an, dass wir eine Eigenlösung Z, des zweiparametrigen Eigenwert- 
problems (4.1) gefunden haben, die zu den Eigenwerten A und 4% gehört. Wir 


546 A. Rubinowicz 


können dann f. auch als eine Lösung eines einparametrigen Eigenwertproblems (1.1) 
: ip 


ansehen, in dem entweder 
(a) q (x) = qo (x) — FE 03 (x), 0 (x) = Qı (x), A= AQ? oder 
(b) q (x) = مو‎ ) — AP o, (9), e (z) = es (9), 4 = 9 


ist. Vorausgesetzt wird dabei, dass in den beiden einparametrigen Eigenwertproble- 
men (a) und (b) das gleiche Grundgebiet wie in dem zweiparametrigen Eigenwert- 
problem (4.1) vorgegeben ist und die gleichen Randbedingungen (1.2a), (1.2b) vor- 
geschrieben werden’). Richtig ist aber auch die Behauptung, dass jede Lösung fe 
des einparametrigen Eigenwertproblems (a) oder (b) zugleich auch eine Lösung des 
zweiparametrigen (4.1) und des einparametrigen (b) bzw. (a) ist. 

Will man also feststellen, für welche Eigenwerte 4% und 2® das zweiparametrige 
Eigenwertproblem (4.1) lósbar ist, so genügt es bei beliebig vorgegebenem AO bzw. 
AO zu ermitteln, für welche Eigenwerte A = A bzw. A = 4? das Eigenwertproblem (a) 
bzw. (b) lósbar ist. Uns interessiert insbesondere der Fall, wo das Eigenwertproblem 
(a) bzw. (b) nicht entartet ist und ein diskretes Eigenwertspektrum besitzt, falls 2 
bzw. AU? in gewissen Bereichen beliebig vorgegeben werden. 

Im Falle des Eigenwertproblems (a) gehórt dann z. B. zu jedem Wert des Para- 
meters AO? eine diskrete Folge von Eigenwerten A = AU. In der A, AO-Ebene 
werden somit die zusammengehörigen Werte der Eigenwertparameter AY und 4% 
durch „Eigenwertkurven“ (Collatz 1949) gegeben, die alle zusammen eine Kurven- 
schar bilden. Eine solche Kurvenschar legt aber im allgemeinen die zusammengehó- 
rigen Werte der beiden Eigenwertparameter nicht nur im Falle des einparametrigen 
Problems (a) fest, sondern auch im Falle des einparametrigen Eigenwertproblems (b) 
oder auch des zweiparametrigen (4.1). Aus der in Rede stehenden graphischen : Dar- 
stellung der zusammengehörigen Eigenwerte von AU und 42) folgt insbesondere, 
dass zugleich mit dem Eigenwertproblem (a) auch das Eigenwertproblem (b) im 
allgemeinen diskrete 20 Eigenwerte bei beliebig vorgegebenen A(? aufweisen wird. 

Der Fall dass die Eigenwertkurven, die sich beim Lösen des einparametrigen 
Eigenwertproblems (a) für verschiedene A(?-Werte ergeben, nicht unbedingt ein 
diskretes AD. -Eigenwertspektrum zur Folge haben, tritt nur dann ein, wenn die Eigen- 
wertkurven in der 4%, A®.Ebene durch Parallele zur A®-Achse gegeben sind. In 
diesem Falle erhalten wir stets die gleichen 4% -Eigenwerte unabhängig von den Wer- 
ten AO = const. 

Die Schar der Eigenwertkurven wird im allgemeinen nur in einem Teilgebiet 
der 29) A®.Ebene definiert sein, je nachdem für welche Werte von A und 4? 
Eigenlósungen des ur agen Eigenwertproblems (4.1) vorhanden sind. 


3) Auch aus der. Tatsache, dass für reelle A-bzw. A()-Werte das einparametrige Eigenwert- 
problem (a) bzw. (b) als solches nur reelle Eigenwerte A) bzw. A(2) besizt, kann man schliessen, dass _ 


sobald einer der beiden Eigenwerte AC) oder A9) beim zweiparametrigen Eigenwertproblem (4.1) reell 
st, auch der andere reell sein muss. 
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Wir wollen annehmen, dass die einzelnen Kurven der Eigenwertkurvenschar 
durch einen ganzzahligen Index numeriert sind. 

Man kann sich nun fragen: 

(æ) In welchen Sinne kann man zwei zu verschiedenen zusammengehörigen Eigen- 
wertpaaren Al, A® gehörige Eigenfunktionen als orthogonal bezeichnen? 

(8) In welcher Weise kann man aus der Gesamtheit der Eigenlösungen des zwei- 
parametrigen Eigenwertproblems (4.1) vollständige Systeme von Eigenfunktionen 
auswählen, die im Sinne der üblichen 'Normierungs- und Orthogonalitätsbedingung 
(1.2) der linearen einparametrigen Eigenwertprobleme orthonormal sind? 

Um zu den Orthogonalitätsbedingungen für die Lösungen des zweiparametrigen 
Eigenwertproblems (s.1) zu gelangen und die Frage (x) zu beantworten, multiplizieren 
wir die Gleichung (4.1) für f = f, mit f, und subtrahieren davon den mit f, multipli- 
zierten konjugiert komplexen Wert der Differentialgleichung (4.1) für .كر‎ Unter der 
Voraussetzung, dass die in Betracht kommenden Eigenwerte reell sind, erhalten wir so 


Ta G= = — A — A) e (z) + An? — Am) os Q0] ffs (4.4) 
wo formell ebenso wie in (1.2b) G,, durch 
df, م‎ dfm 
wn (en) 
definiert ist. Integration der Beziehung (4.4) über das Grundgebiet (x,, xa) ergibt dann 


PFA 169 — 39) o. (9 + AP — 19) s, (4.5) 


x 


vorausgesetzt, dass die Randbedingung (1.2a) 
lim Gom = lim Gin 


x>X, xX, 
erfiillt ist, so dass das Integral (4.5) konvergiert. 

- Analog wie wir es in (3.5) getan haben, können wir auch hier zwei Eigenfunktionen 
f, und f,, des zweiparametrigen Eigenwertproblems (4.1) als zueiriander orthogonal 
bezeichnen, falls sie die Bedingungen (4.5) erfiillen. Ganz ebenso wie im Falle des 
quadratischen Eigenwertproblems verhindert aber auch hier das Auftreten der Eigen- 
werte in der Orthogonalititsbedingung (4.5) eine einfache Berechnung der Ent- 
wicklungskoeffizienten 'c, der Entwicklung 


f (z) = 25 Sa @) e, 


einer Funktion f(x), selbst wenn eine solche Entwicklung tatsächlich vorhanden ist. 
Um die Beantwortung der Frage )6( vorzubereiten, wollen wir zeigen, dass man 
sehr leicht von zweiparametrigen Eigenwertproblemen zu einparametrigen zurück- 


١ gelangen kann, für die die üblichen Normierungs- und Orthogonalitätsbedingungen 
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(1.2) verwendbar sind. Spezielle Fälle von einparametrigen Eigenwertproblemen, die 
in zweiparametrigen enthalten sind, sind die Eigenwertprobleme (a) und (b). Um 
allgemeinere Fälle zu erhalten, in denen die Orthogonalitätsbedingung (4.5) die Gestalt 
der Bedingung (1.2) annimmt, muss man voraussetzen, dass zwischen zwei beliebigen 
Eigenwertpaaren AD, 2 und 4%, A® eine Beziehung von der Gestalt 


A (AD — 4) = B (a — 40) (4.6) 


besteht. In der Orthogonalitätsbedingung (1.2) wird dann gemäss (4.5) die Gewichts- 
funktion o(x) bis auf eine multiplikative Konstante durch 


e (x) = Be, (x) + Aes (x) (4.7) 
gegeben. Die Orthogonalitätsbedingung (4.5) geht dann über in 


f L A, Bor (a) + Aos 69] dx = 0. (48) 


Spezielle Fälle von (4.6) sind A = 0 d.h. A? = const sowie B = 0 d.h. A = const. 
Auf Grund von (4.7) wird dann o (x) = Bo, (x) bzw. o (x) = A وم‎ (x), so dass das zwei- 
parametrige Eigenwertproblem (4.1) in das einparametrige Eigenwertproblem (a) 
bzw. (b) übergeht. 

Im allgemeinen Falle bedeutet jedoch (4.6) dass 


AMM — Bà? — C (A,B, C = const) (4.9) 


ist, also eine lineare Beziehung zwischen den Eigenwertparametern AU und 22 
besteht. In der A, A®-Ebene entspricht der Beziehung (4.9) eine Gerade. Ihre 
Schnittpunkte mit der Schar der Eigenwertkurven legen eine Folge von diskreten 
Werten der beiden Eigenwertparameter A und 4? fest. Die Eigenfunktionen des 
zweiparametrigen Eigenwertproblems (4.1), die zu diesen Paaren von Eigenwert- 
parametern gehóren, sind nicht nur im Sinne der Beziehung (1.2) zueinander ortho- 
gonal, sondern kónnen auch im Sinne dieser Beziehung-normiert werden. Eine natür- 
liche Numerierung der zusammengehórigen Paare von Eigenwertparametern A? 
und A? ergeben die Nummern der Eigenwertkurven auf denen sie liegen. 

Die Folgen dieser zusammengehórigen Eigenwertpaate A und A® und die 
entsprechenden Eigenfunktionen kann man auch aus einparametrigen Eigenwert- 
problemen ermitteln, wenn man i.. dem zweiparametrigen Eigenwertproblem (4.1) 
den Parameter A bzw. 4? mit Hilfe der Beziehung (4.9) eliminiert. Man erhält 
auf diese Weise die nachstehenden beiden einparametrigen Eigenwertprobleme 


d d 
E - ac هد - وميه ع‎ (ai ZONE 
bzw. 


d d, G | 
ERA (x) = sF E (x) — et (x) — 49 (e (x) + = 01 هدم | زه‎ 


^ 
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Sie sind einander vollständig äquivalent, da sie die gleichen Eigenfunktionen und 
auf Grund von (4.9) auch die gleichen Paare von zusammengehörigen Eigenwerten A 
und A? ergeben. | 

Die Wahl der in (4.9) auftretenden Konstanten A, B, C wird man Einschränkun- 
gen unterwerfen, jenach den Forderungen, die man an das System der Eigenfunktionen 
stellt, das sich mit Hilfe von (4.9) ergibt. Im allgemeinen wird man z.B. die Schnitt- 
gerade (4.9) so legen, dass das durch (4.7) gegebene o(x) im Grundbereiche positiv ist. 
um eine Normierung der Eigenfunktionen auf + 1 zu ermöglichen. Eine weitere 
Einschränkung in der Wahl der Koeffizienten A, B, C wird im allgemeinen die Forde- 
rung )6( ergeben, die die Vollständigkeit des Systems der Eigenfunktionen verlangt. 
Ein solches System von orthonormalen Eigenfunktionen muss z.B. eine unendliche 
Anzahl solcher Funktionen enthalten. Die durch (4.9) gegebene Gerade muss daher, 
falls die Eigenwerte nicht entartet sind, unendlich viele Eigenwertkurven schneiden, 
da ja jedem Schnittpunkte ein zusammengehöriges Paar der Eigenwerte 2, A des 
Eigenwertproblems (4.1) entspricht. Ob man in der Tat beim Erfülltsein dieser not- 
wendigen Bedingung ein vollständiges System von orthonormalen Eigenfunktionen 
erhält, muss man in jedem Spezialfalle besonders untersuchen. 

Soweit mir bekannt ist, lassen sich Beispielle von linearen zweiparametrigen 
Eigenwertproblemen, die für die Quantenphysik von Bedeutung sind, nicht anführen‘). 
Mathematische Beispiele für solche Probleme (4.1) kann man jedoch in allen Fällen 
aus einparametrigen Eigenwertproblemen (1.1) erhalten, sobald sie sich ‚wie in $2 
angegeben, umordnen lassen und dabei lineare umgeordnete Eigenwertprobleme 
ergeben. 

Dies gilt insbesondere für die mit Hilfe der Polynommethode lösbaren Eigen- 
wertprobleme, die in lineare Eigenwertprobleme umgeordnet werden können?). 
Der Mechanismus, der es bewirkt dass wir für diese spezielle Klasse von Eigenwert- 
problemen im allgemeinen Scharen von Eigenwertkurven erhalten, ist sehr einfach. 
Um Eigenwerte zu bekommen muss man, wie bereits in 8 1 erwähnt wurde, im 
Falle eines diskreten Eigenwertspektrums in den in den Lösungen auftretenden hyper- 
geometrischen Funktionen >F; (a, b; c;£) bzw. F, (a;c; £) einen von den beiden 
Parametern a oder b bzw. den Parameter a gleich einer negativen ganzen Zahl setzen. 
Eine solche Festsetzung der Eigenwerte kann selbstverständlich nur dann erfolgen, 
wenn die genannten Parameter a oder b von dem Eigenwertparameter abhängen. Soll 
daher eine Beziehung zwischen und 4% und 4% bestehen, so muss der Parameter a 
oder b, der einer negativen ganzen Zahl gleichgesetzt wird, gleichzeitig eine Funktion 
von A® und A® sein. Bei den mit Hilfe der Polynommethode lösbaren Eigenwert- 


men ا‎ Fa do E ري‎ BE ص‎ E S 
4) Zweiparametrige Eigenwertprobleme die in der Technik auftreten führt Collatz (1949) an. 
5) Lassen sich lineare einparametrige Eigenwertprobleme in quadratische umordnen (vgl. $3), 
so ergeben sie zweiparametrige Eigenwertprobleme, die in dem einen Eigenwertparameter linear, in dem. 
anderen jedoch quadratisch sind. Mit solchen zweiparametrigen Eigenwertproblemen wollen wir uns 
hier nicht befassen, wenn sich; auch Beispiele für solche Probleme angeben lassen, die mit Hilfe der 


Polynommethode lösbar sind. 
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problemen werden somit die Eigenwertkurven durch eine Beziehung zwischen den 
beiden Eigenwertparametern A® und 2® und der Polynomquantenzahl gegeben. 
Da die Eigenwertparameter Funktionen der entsprechenden Eigenwertquantenzahlen 
sind, müssen diese Beziehungen äquivalent sein den Relationen, die zwischen den 
Eigenwertquantenzahlen und der Polynomquantenzahl bestehen. Jedem Wert der 
Polynomquantenzahl entspricht auf diese Weise eine Eigenwertkurve. Die Polynom- 
quantenzahlen numerieren sozusagen ganz von selbst die Eigenwertkurven. 

Im Falle des Eigenwertproblems der zugeordneten Kugelfunktionen (1.3) er- 
geben sich die Gleichungen der Eigenwertkurven aus der Beziehung 1 = m + n (1.9), 
wenn man hier gemäss (1.8) und (2.2) die Eigenwertquantenzahlen / und m durch 
die entsprechenden Eigenwerte 4% = / (l + 1) und AO = — m? ausdrückt. Man 
erhält so, wenn J und m als positiv vorausgesetzt werden, für die Schar der Eigen- 
wertkurven die Darstellung 


— 1 + V20 TS == + n. (4.10) 


Im Falle des Eigenwertproblems (1.4) der Radialfunktionen eines Einelektro- 
nenatoms werden die Eigenwertkurven aus der Relation n = n, + I+ 1 (1.11) 
erhalten, wenn man berücksichtigt, dass gemäss (1.12) und (2.4) die Eigenwert- 
parameter durch A® = — 1/n2a? und 4% = — l1 (l + 1) gegeben werden. Sie 
werden dargestellt durch 


V — 139 ad =n, ++ V— 19 + 1. (4.11) 


Den Punkten der durch die Polynomquantenzahl n gekennzeichneten Eigen- 
wertkurven (4.10) entsprechen dabei die zugeordneten Kugelfunktionen P(x), bei 
denen die Differenz der beiden Indizes |—7» gleich ist der betreffenden Polynom- 
quantenzahl n. 

Den Punkten der mit der Polynomquantenzahl n, numerierten Eigenwertkurve 
(4.11) sind hingegen die Radialfunktionen R,, (x) zugeordnet, für die die Differenz 
der Indizes n—/ durch die um 1 vermehrte Polynomquantenzahl n, d.h. durch n,4-1 
bestimmt wird. 

Wir wollen uns auch an Beispielen überzeugen, wie die Normierungs- und Ortho- 
gonalitätsbedingungen (4.8) mit den durch (4.7) gegebenen Gewichtsfunktionen aus- 
sehen, insbesondere, da sie in der mathematischen Literatur nicht bekannt zu sein 
scheinen. Im Falle des Eigenwertproblems der zugeordneten Kugelfunktionen ist 
o,(x) = 1 (1.5) und اوم‎ = 1/(1—x?) (2.3). Die entsprechenden Normierungs- und 
Orthogonalitatsrelationen (4.8) lauten daher 


+1 
1 PY" (x) PF (x) Ez | dx = Sn! « (4.12) | 
-1 


Dem Kroneckerschen Delta-Symbol haben wir dabei als Indizes die De 
zahlen der Eigenwertkurven beigefügt, auf denen die den beiden auftretenden zuge- 
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ordneten Kugelfunktionen entsprechenden Eigenwertpaare AD, 4% liegen. Das 
Integral (4.12) muss nämlich verschwinden, falls die Eigenwertpaare AU, 4% der 
beiden Eigenfunktionen auf verschiedenen Eigenwertkurven liegen und muss gleich 
1 sein, falls sie in einem Punkte einer Eigenwertkurve zusammenfallen. 

Um die J, m- bzw. P, m'-Werte festzulegen, für die die Bedingung (4.12) erfüllt 
ist, beachte man, dass l = m+n (1.9) ist, und daher mit Rücksicht auf A = L (1+1) 
(1.8), 4% = —m? (2.2) sich aus (4.9) die Beziehung 


A41 (1l + 1) + Bm? = C (4.13) 


ergibt. Wir erhalten daher mit Hilfe von (4.13) im Falle A+ B= 0 für l und m die 
nachstehenden Ausdrücke 


A A? A 11 C+} 
m = — —  — (2 n + 1) + 2 d 5 +A 
JEDE V [aas TM RE. 
Í = m + n. (4.14) 
Im Falle A+B = 0 gilt hingegen 
_ n(n4 1) I G 
لسن‎ ne 
| — m + n. (4.15) 


Für /', m’ ergeben sich entsprechende Ausdrücke mit L, m na 1 Lemon. 
Aus (4.14) sind leicht die Bedingungen abzulesen, dass das durch die Relationen 

(4.14) ausgewählte System von Eigenfunktionen unendlich viele solche Funktionen 

enthält. Für grosse n-Werte ergibt sich aus (4.14) der nachstehende Ausdruck für m 


A B 
Damit die in (4.16) auftretende‘ Quadratwurzel reell ist, muss 
$ <0 T (4.17) 


sein. Eine weitergehende Einschränkung für den Quotienten B/A ergibt sich aus der 
Forderung, dass m;20 sein muss. Der unseren Überlegungen zu Grunde liegende 
Ausdruck (1.7) für die zugeordneten Kugelfunktionen wurde nämlich unter der 
Voraussetzung m 7z 0 und 1720 abgeleitet. Die bei der Berechnung von (4.14) benutzte 
Relation (1.9) gilt auch nur unter dieser Voraussetzung. 
In dem in (4.17) enthaltenen Spezialfalle ” 
B 


; -1<7<0 (4.18) 
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ist jedoch die in (4.16) auftretende Quadratwurzel kleiner als 1 und der vor der Klammer 
in (4.16) stehende Ausdruck A/2(4+B) positiv. Im Falle (4. 18) ergibt sich daher 
aus (4.16) für m ein negativer, also unbrauchbarer Wert. | 

In dem noch zur Untersuchung verbleibendem Spezialfalle von (4.17) 


I (4.19) 


ie nd 

hat hingegen die Quadratwurzel in (4.16) einen Absolutbetrag, der grósser als 1 ist, 
während der hier vor der Klammer auftretende Ausdruck A/2(4+B) negativ ist. 
Wir müssen daher bei der in (4.16) auftretenden Quadratwurzel das negative Vor- 
zeichen berücksichtigen, um für m einen positiven Wert zu erhalten. Aus (4.19) 


folgt jedoch 


TS = < 0. (4.20) 


Schreibt man die Gleichung (4.9) für die Schnittgerade in der Form 


A C 
2 T A — B (4.21) 
so erkennt man, dass (4.20) fordert, dass die Schnittgerade in der A, A®.Ebene 
mit der AU-Achse einen Winkel zwischen 0 und — iz einschliesst. 

Dieses Ergebnis ist verständlich, wenn wir bemerken, dass gemäss (1.8) AD = 
= l(l--1) einen positiven und gemäss (2.2) 4% = — m? einen negativen Wert hat 
und daher die Eigenwertkurven in dem vierten Quadranten der 49) 4@-Ebene - 
verlaufen. Sie werden hier mit Rücksicht auf 7 = m+n (1.9) durch die Gleichung 


E A® (RFE EVADE) (4.22) 


gegeben. Für grosse A®.Werte gehen sämtliche Eigenwertkurven, wie aus (4.22) 
sich ergibt und auch aus der Fig. 1 zu entnehmen ist, asymptotisch in Gerade über, 
die zur gestrichelten Winkelhalbierenden der positiven A®-und der negativen 
AO-Achse parallel sind. | 

Der durch (4.15) wiedergegebene Fall A+B — 0 entspricht, wie aus (4.21) 
ersichtlich ist, dem Grenzfall, wo der Neigungswinkel der Schnittgeraden gerade — 1z 
beträgt. 

Aus (4.14) kann man auch die: Bedingung entnehmen, die erfüllt sein muss, 
damit die Schnittgerade (4.21) alle Eigenwertkurven schneidet. In diesem Falle muss 


(4.14) auch für. n = 0 dh. 


2 xe] ) . (4.23) 


n= a y 
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reell und positiv sein. Damit die in (4.23) auftretende Quadratwurzel reell ist, muss 


A? + AC(A +B)>0 d.h. 


A?JB? 0 


En Da die linke Seite von (4.24) gemäss (4.20) positiv ist, bedeutet (4.24) dass der 
Schnittpunkt — C/B der Geraden (4.21) mit der 2®-Achse je nach dem Werte von 


0 10 20 30 40 50 60 
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2 2 2 
! ` 
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Fig. 1. Eigenwertkurven des zweiparametrigen Eigenwertproblems der zugeordneten Kugelfunktionen. 
Die strichlierte Gerade entspricht der Winkelhalbierenden zwischen der positiven A@)- und der nega- 
tiven A(?)-Achse. j 


A/B, der der Einschränkung (4.20) unterliegt, nicht zu hoch auf der positiven 4®-Achse 
liegen darf. Aus (4.23) entnimmt man ferner die Tatsache, dass mit Rücksicht auf 
(4.19) der Ausdruck unter dem Wurzelzeichen für C/A>0 kleiner ist als 1 und wir 
daher für m zwei verschiedene Werte erhalten. Für C/A <0 ist jedoch dieser Ausdruck 
grösser als 1, so dass nur ein m-Wert zulässig ist. Mit Rücksicht darauf, dass C/A 
den Schnittpunkt der Geraden (4.9) mit der A? -Achse bezeichnet, bedeutet dies. 
dass sobald die Schnittgerade (4.9) die negative A -Achse schneidet, sie mit der durch 

= 0 gegebenen Eigenwertkurve nur einen Schnittpunkt gemeinsam hat, während 
sich im Falle, wo sie die positive A-Achse schneidet, zwei Schnittpunkte ergeben. 
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Im Falle der Radialfunktionen R,, (x) ist o, (x) = x? (1.4) und o, (x) = 1 (2.5). 
Daher wird die Normierungs- und Orthogonalitätsbedingung (4.8) gegeben durch 


f Ry, (2) Ryp (x) [Bx? + A] dx = ónn;. 
0 


da die, die Eigenwertkurven numerierende Polynomquantenzahl jetzt mit n, bezeichnet 
wird. 

Um die Eigenwertquantenzahlen n,| in ihrer Abhängigkeit von der Polynom- 
quantenzahl n, zu ermitteln, bemerken wir, dass n = n, +I + 1 (1.11) ist und daher 
mit Rücksicht auf A? = — 1/n?a2 (1.12), A9 = — l (l + 1) (2.4), sowie (4.9) die 
Eigenwertquantenzahlen n und / in ihrer Abhängigkeit von n, aus den Gleichungen 


- وك كد‎ + Bl + - 6 E 


n=n, titl 


zu berechnen sind. Auf diese Weise ergeben sich für n und / algebraische Gleichungen 
vom vierten Grade. Analoge Gleichungen gelten für n',/',n,'. Es mag daher auf 
die weitere Diskussion dieses Falles verzichtet werden. 

Es sei noch darauf hingewiesen, dass sich in ähnlicher Weise auch dreiparametrige 
Eigenwertprobleme behandeln lassen. Ein Beispiel für ein solches kann aus (2.8) 
erhalten werden. Setzt man hier AP = l (I + 1), 4% — — 1 m?, A9 = — 1 n?, so 
erhält man das dreiparametrige Eigenwertproblem 


d df 1 1 
et “yay ال‎ ae Ck) ee S10) ==) (3) = 
Ar | 1 — 92) + Ar م‎ 0. 


Die Verallgemeinerung der bei den zweiparametrigen Eigenwertproblemen durch- 
geführten Überlegungen auf den dreiparametrigen Fall bietet keine Schwierigkeiten. 
Auch wenn die Differentialgleichung des Eigenwertproblems selbstadjungiert ist, 
müssen selbstverstándlich die Eigenwerte nicht unbedingt reell sein. Ferner treten 
an Stelle der Eigenwertkurven in der A", A®.Ebene hier in einem dreidimensio- 
nalen AP, 2, A®-Raume- Eigenwertflächen auf, die den verschiedenen Polynom- 
quantenzahlen entsprechen. Man erhält zweiparametrige Eigenwertprobleme aus 
dreiparametrigen, falls man in dem 4%, 22, 2-Raume die Flächenschar der Eigen- 
wertflächen durch eine Ebene schneidet. Wird die Schar der Eigenwertflächen durch 
eine Gerade geschnitten, so ergeben sich Sätze von Eigenfunktionen, die im Sinne der 
üblichen Orthogonalitätsbedingung (1.2) zueinander orthogonal sind und die auch 
als Lösungen einparametriger Eigenwertprobleme aufgefasst werden können. 

Aus den mit Hilfe der Polynommethode.lósbaren Eigenwertproblemen kann 
man höchstens dreir arametrige Eigenwertprobleme erhalten, da j ja aus einem solchen 
ursprünglichen Eigenwertproblem laut 82 nur höchstens zwei umgeordnete Eigen- 
wertprobleme herstellbar sind. 
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85. Beziehungen zur Faktorisierungsmethode 


Die Eigenfunktionen, die mit Hilfe der Faktorisierungs methode (Schródinger 
1940, 1941a und b, Infeld 1941, 1947, Infeld-Hull 1951) erhalten werden, hángen in 
der Regel®) ausser von den Eigenwerten A, des Eigenwertparameters auch noch von 
einem Parameter m ab, der gewóhnlich als ganzzahlig angenommen wird. Wir kónnen 
sie daher mit 


u = u (Àp: m, x) (5.1) 


bezeichnen. Dabei wird ohne nähere Begründung vorausgesetzt, dass die Eigenwerte An 
der Eigenfunktionen (5.1) unabhängig vom Zahlenwert des ganzzahligen Parameters m 
stets einer Teilfolge der gleichen Folge von Zahlen zu entnehmen sind. 


Der in der Faktorisierungsmethode auftretende Parameter m ist vielfach physikali- 
schen Ursprungs. In der Regel handelt es sich nämlich bei faktorisierbaren Eigen- 
wertproblemen um Probleme, die durch Separation von Variablen aus einer partiellen’ 
Differentialgleichung gewonnen werden und m bedeutet dann gewóhnlich eine Eigen- 
wertquantenzahl eines Eigenwertproblems, das vor Inangriffnahme desjenigen Pro- 
blems gelöst werden muss, das m als Parameter enthält. 

Vom physikalischen Standpunkt ist, wenigstens in vielen solchen Fällen, die 
Unabhängigkeit der A,-Werte vom Parameter m plausibel, insbesondere wenn, wie 
oben angedeutet, der Parameter m physikalischen Ursprungs ist. Beim Eigenwertpro- 
blem (1.3) der zugeordneten Kugelfunktionen bestimmen die Eigenwerte A = /(/4- 1) 
(1.8), wenn man sie mit dem Faktor A? versieht, die Quadrate des gesamten Bahn- 
impulsmomentes. Die Tatsache, dass ohne Rücksicht auf den Wert von m die Eigen- 
werte stets ein und derselben Folge von Zahlen zu entnehmen sind, ist physikalisch 
plausibel,da sie ja nichts anderes bedeutet, als dass bei gegebener Projektion mh des 
Bahnimpulsmomentes das Quadrat des gesamten Bahnimpulsmomentes durch jeden 
Wert I(l--1)h? gegeben sein kann, falls nur 7 der Bedingung |m|</ entspricht: 

Im Falle des Eigenwertproblems (1.4) der Radialfunktionen eines Einelektronen- 
atoms bestimmen die Eigenwerte A (1.12) die Energieniveaus dieses Atoms. Der 
Parameter l legt hingegen das gesamte Impulsmoment der Bahn fest. Dass die Folge 
der Eigenwerte unabhängig vom Werte von / stets der gleichen Zahlenfolge (1.12) 
angehórt, ist physikalisch plausibel, wenn man beachtet, dass die Energie der Bahn 
nur durch die ‚on J unabhängige grosse Achse nicht aber durch die durch / mitbes- 
timmte Exzentrizität der Bahn festgelegt wird. | | 

Bisher wurde jedoch noch nicht versucht vom mathematischen Standpunkt 
den oben angeführten sehr frappanten Tatbestand zu erklären. Um zu zeigen in welcher 
Richtung eine solche mathematische Erklärung zu suchen ist, bemerken wir, dass 


EEE en a gg ب تس سس‎ NE Ta BT eee ee E COO IO SS TUI 


| 6) Eine Ausnahme bildet z. B. das Eigenwertproblem des linearen harmonischen Oszillators. 
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wenigstens alle im Infeldschen Sinne faktorisierbaren Eigenwertprobleme’) (Króli- 
kowska 1958, Rubinowicz M) mit Hilfe der Polynommethode gelóst werden kónnen. 
Für die weitere Klasse der im Schródingerschen Sinne faktorisierbaren Eigenwert- 
probleme wurde dies zwar nicht bewiesen, es ist jedoch kein in diesem Sinne faktorisier- 
bares Eigenwertproblem bekannt, das nicht mit Hilfe der Polynommethode gelöst 
werden kónnte. Wir wollen daher bei unseren weiteren Überlegungen nur Eigen- 
wertprobleme in Betracht ziehen, die mit Hilfe der Polynommethode lósbar sind. 

Wir behaupten nun dass der eingangs dieses Paragraphen erwühnte Tatbestand 
erklürlich ist, sobald wir den Parameter m in den faktorisierbaren Eigenwertproblemen 
als die Eigenwertquantenzahl eines aus dem gegebenen Eigenwertproblem umgeordne- 
ten Eigenwertproblems ansehen und dabei voraussetzen, dass faktorisierbare Eigen- 
wertprobleme auch mit Hilfe der Polynommethode behandelt werden kónnen. 

Die Richtigkeit dieser Behauptung können wir zunächst an den beiden speziellen 
Eigenwertproblemen feststellen, deren wir uns zur Erläuterung der allgemeinen Sätze 
bedient haben. Im Falle des Eigenwertproblems (1.3) der zugeordneten Kugelfunktio- 
nen wurde die Eigenwertquantenzahl.des umgeordneten Problems mit m und im 
Falle des Eigenwertproblems (1.4) der Radialfunktionen mit 7 bezeichnet. Wie der 
Vergleich mit Schrödinger (1940, 194la und b) und Infeld-Hull (1951) zeigt 
spielen diese Zahlen in der Tat die Rolle des oben mit m bezeichneten Parameters. 
Aus der Faktorisierungstabelle bei Infeld-Hull (1951, S. 66) kann man aber auch 
entnehmen, dass in allen im Infeldschen Sinne faktorisierbaren Eigenwertproblemen m 
als die Eigenwertquantenzahl eines entsprechenden umgeordneten Eigenwertproblems 
angesehen werden kann. 

Die Unabhängigkeit der Eigenwerte A von dem Parameter m folgt in unseren 
beiden Beispielen aus den Beziehungen, die zwischen der Polynomquantenzahl und 
‚den Eigenwertquantenzahlen des ursprünglichen und des umgeordneten Eigen- 
wertproblems bestehen. 

Im Falle des Eigenwertproblems (1.3) der zugeordneten Kugelfunktionen lautet 
diese Beziehung l = m + n (1.9), wobei 2 = I(l + 1) (1.8) d.h. 


à = (m+n)(m+n+1), n=0,1,2,... (5.2) 


ist. Ist nun m eine ganze Zahl, so folgt aus (5.2) unmittelbar dass die A-Werte ohne 
Rücksicht auf den Wert von m stets der Folge der Zahlen A = 0, 2, 6, 12, ... zu 


entnehmen sind. Unsere Behauptung ist aber mutatis mutandis auch richtig, wenn m 
die allgemeinere Gestalt 


m = m'+Am, Am = const, m' = 0,1,2, ... 
hat, also der Parameter m sich nur um ganze Zahlen ändern kann. Die A-Werte ge- 
Se i ae m 
?) Bezüglich der Definition der im Schródingerschen und im Infeldschen Sinne faktorisienbaren 
Eigenwertprobleme vgl. Rubinowicz M. 
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hören dann der Zahlenfolge 


À = Am(l + Am), (1 + Am)(2 + Am), (2 + Am)(3 + Am), ... 
an. 
Im Falle des Eigenwertproblems der Radialfunktionen gilt n = n,4-l--1 (1.11), 
wobei A durch (1.12) gegeben wird. Es ist demnach 


1 
an, + 1+ 1 


Analog wie im Falle der zugeordneten Kugelfunktionen schliesst man auch hier, mit 
Rücksicht auf die Tatsache, dass die Polynomquantenzahl n, ganzzahlig ist, dass 2 
ohne Rücksicht auf den Wert von / stets eine Teilfolge der gleichen Zahlenfolge ist. 
Dies gilt nicht nur wenn J, wie dies die Quantentheorie fordert, durch eine ganze 
Zahl gegeben wird, sondern auch in dem Falle, wo Í in der Gestalt I = 7 + AL (Al = 
= const, /' = 0, 1, 2, ...) vorausgesetzt wird. 


A= 


Der Nachweis, den wir an den beiden spezielien Beispielen durchgefiihrt haben, 
dass die Folge der Eigenwerte A unabhängig vom Werte des ganzzahligen Parameters m 
ist (oder von der ganzzahligen Differenz der m-Werte, im Falle nicht ganzzahliger m) 
zeigt die Richtung in der der Beweis im allgemeinen Falle zu suchen ist. Es geniigt 
nachzuweisen, dass bei den mit Hilfe der Polynommethode lösbaren Eigenwert- 
_problemen die Eigenwertquantenzahl 7 des ursprünglichen Eigenwertproblems bis 
auf eine additive Konstante durch die Summe oder Differenz der Polynomquantenzahl n 
und der Eigenwertquantenzahl m des umgeordneten Eigenwertproblems, also in der 
Gestalt 
t = +n+m-const (5.3) 


darstellbar ist. Die Eigenwerte des ursprünglichen Eigenwertproblems sind dann durch 
2 = Al) = A( + n + m + const) 


gegeben. Eine Änderung von m um eine ganze Zahl kann dann durch eine entspre- 
chende Änderung der ganzzahligen Polynomquantenzahl n kompensiert werden, 
vorausgesetzt dass n nach dieser Änderung noch immer positiv bleibt. 

Dazu muss noch bemerkt werden, dass es für die Richtigkeit unserer Behauptung 
durchaus nicht erforderlich ist, dass das umgeordnete Eigenwertproblem linear ist. 
Das kommt ja in der Beziehung (5.3) in keiner Weise zur Geltung. Das umge- 
ordnete Eigenwertproblem kann daher auch quadratisch sein, wie dies bei einigen 
Eigenwertproblemen tatsüchlich der Fall ist, die den Infeldschen Typen A und C 
(vgl. Infeld-Hull 1951, S. 66 und 67) angehören. Der Unterschied zwischen den 
Eigenwertproblemen, die nach einer Umordnung lineare bzw. quadratische Eigen- 
wertprobleme ergeben, besteht nur darin, dass bei den ersteren die umgeordneten 
Eigenwertprobleme im allgemeinen im Infeldschen Sinne faktorisierbar sind, während 
bei den letzteren dies ausgeschlossen ist. Die Infeldsche Faktorisierungsmethode 
ist ja doch nur auf lineare Eigenwertprobleme anwendbar. 
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In this paper a charge singlet for the recently discovered D-particle is proposed. 
This particle is included into the scheme of Rayski after a slight modification. Other 
possibilities in the framework of this scheme are also discussed. 


I. Introduction 


The existence of the D-particles may be regarded as something rather established!). 
Since their isotopic properties, strangeness etc. are still unknown, it seemed worth- 
while to make some speculations thereupon. The discussion from the point of view 
of the schemes of Nishijima-Gell-Mann aud of Salam and Polkinghorne has been 
performed recently by Bialkowski and Jurewicz [1]. In the present paper some specu- 
lations are being made in the framework of the schemes of Rayski and Sakata. 


II. The systematics of Rayski 


Since it seems quite sure that the D-particle is a boson, it may be investigated 
whether it can be included into the table of bosons which was given by Rayski in his 
last paper [2]. The basis of his systematization was the formula: 


q— n Tg + fs i Q) 


where n stands for the baryon number, but it is also equal to 1 for the K-mesons. 
If we take n — 1 for the further mesons (with the exception of the pion), the ordering 
is similar to that of the fermions (a neutral singlet,a charge triplet and a charge doublet). 
For the D-particle it is the simplest to assume the charge triplet which in this system- 
atization scheme should give lower values of masses than the doublet. 

It is thus possible that for the quantum numbers Tj, = —1 and Ig = —1, 
0, 1 we will have the three paritcles Dt, D° and D- (in the order of raising masses); 
pri coge, "Under Vat sop RN TS; ME WE UE 

1) T. Yamanouchi, preprint. 
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the corresponding antiparticles D+, D? and D` (different from the particles) may be 
obtained by taking Ts, = 1. For Tj, = —3/2 the only particles would be D° and D^, 
whereas the positively charged D-meson would be the antiparticle to D". 


III. Tests for the isotopic structure in the Rayski scheme 


In order to distinguish between the possibilities of triplet and doublet it is impor- 
tant to study the production and decay of the D-mesons. Since the new quantity مج‎ 
is now taken as a physically meaningful substitute for the phenomenological concept 
of strangeness, we will use Rayski's selection rules (which come from a modification 
of Gell-Mann's rules): AI; = 0 for strong and electromagnetic interactions; A/;, = 
= +) for weak interactions, while AJ;, = + lis for very weak or none interactions. 
Let us take these selection rules as the basis for our considerations on the two possi- 
bilities. 


a) The triplet case: 


The D-mesons may be produced in pion-nucleon and nucleon-nucleon reactions. 
They are allowed to appear in pairs, e. g. 


D---D- +n 
D+ D + n 
D- 4- D? - p (2) 
D+ D* 4- p 


gz D> 


and in some other processes involving at least three particles in final state like the 
reaction i 

D + Z- + Kt 

D- + 29 + Kt 

D+ Z* + K? 

Dt + 2- + K? (8) 


A simple associated production of D-particles and hyperons is, however, forbidden 
in such reactions. The possibleKz-decays are: 


ms dp 


K? + سج‎ (4.1) 

D- 
a | K-+m (4.2) 
Kt + m (A.3) 

p+ 
a | Kat (4.4) 
Kan (4.5) 
Dt). K^ dons | (4.6) 


xt (4.7)‏ ل 
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Rayski's selection rules for strong resp. weak interactions are, of course, fulfilled 
for the production and decay modes above. 


b) The doublet case: 


In this:case production channels for D-mesons are opened. Hyperon-D-meson 
pairs can be now produced in pion-nucleon and nucleon-nucleon reactions. The 
hyperon associated with D-meson can be X or A, but not £. These are some of the 
simplest production modes: 


! pues 
uj wey rud Do + Xo (5) 
D? + 4° 
D-+Z04p 
ptp> |D-+4%+p (6) 
DUNS. 


What regards the Kz-decay modes of the D-particles, there is following distinction 
with respect to the triplet case: The process (4.1) must be replaced by: 


D-— K? + a7 (4.1) 
and the decay (4.5) is now impossible. The remaining decay modes ra with 


(4.1’) fulfil Rayski’s selection rule for strong interactions. Since this is very strange, 
a comparison with experiment in necessary. 


c) Comparison with experiment 


There are three reported observations which are now being interpreted as involv- 
ing a strong interaction of the D-particles with nuclei. The reaction occurring is'as- 
sumed to be 


D- + 2p > Z+ ل‎ 4° t (7 
for the case reported in [3], and 
D- + n— K + A (8) 


for the two events reported in the papers [4] and [5]. It may be discussed whether 
the selection rule 4/,, = 0 is fulfilled in these both reactions (which are regarded 
as “hot“). Since all values of يو[‎ for the particles occurring (except of the D-meson) 
are given in the paper of Rayski [2], from our assumption that the processes (7) and (8) 
are strong follows that we have to take وول‎ = 0 for the D-meson. This is in contra- 
diction with our previous statements. This apparent contradiction can be abandoned 
if slight modifications of the theory presented would be performed. It is thus possible 
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to take n = 2 and T4, = —1 for D-mesons in order to obtain agreement with expe- 
riment. It is, however, difficult to interpret this, since n is usually assumed as baryon 
number. A simple way to extend the validity of Rayski's scheme to the D-particle 
is to adopt the idea of a charge singlet from the Sakata model. 


IV. The Sakata model 


It is well known that in the Sakata model [6] we have two single charged mesons 
for which there are empty spaces in the Gell-Mann scheme. These are the wt and w~ 
mesons, which may be represented as follows: 


ot = 75 (pn — ap) AA 


1 
(DES ya np) AA 
They have isotopic spin zero and the strangeness of the first of them is +2 ( — 2 for 
the second). They may be included into the systematization of Rayski in a straight- 
forward manner if the meaning of the formula (1) is altered for bosons, namely if we 
take always n = 0 for them. In this case we have the following table for the bosons ; 


Table 1 


The ordering, like in Rayski's original paper [2] is such that masses are raising accord- 
ing to their distance from the centre. We see that the w-meson which from now will 
be identified with the D-particle is the heaviest meson among the known ones. The 
existence of further mesons is rather doubtful since in the further multiplets there 
would appear at least one multiple charged elementary particle and this is too vague 
(at least presently) to be assumed. It seems that the requirement that elementary 
particles be single charged is sufficient tó cut off from the a priori infinite table of 
elementary particles only a small part. 

Let us compare our assumption with the experiment. The processes (7) and (8) 
are now strong, since we have Al,, = 0 in both of them. Further we have Al, = 4 
in the Kz-decay modes of the D-particles?) 

Even if new experimental results would contradict the consequences which 
have been drawn here from the reactions (7) and (8), it is quite sure that one will 


2) The D*— meson decays into a Kr-pair, whereas the D-meson decays into the Kz-pair. 
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be able to distinguish between the possibilities of e. .ع‎ a triplet of D-particles in the 
Rayski scheme, and a singlet in the modified Rayski scheme. Associated production 
of D-particles and Z hyperons is now possible, which makes a clear difference between 
the singlet assumed now and the two other assumptions discussed previously. Instead 


of the production modes (5) and (3) we have now only the possibilities: 
E=+ Dt 


= +>] = 
nja eM Ke (9) 


and instead of (6) there is in the singlet case: 
pI pr DEE p-r = (10) 


It seems that the existence of different exit channels in each case gives us an oppor- 
tunity to prove experimentally what one of the presented possibilities for isotopic 
structure of the D- particle is valid. Let us consider e. g. the case of p+p reactions 
which are of practical importance. In these reactions the final products may be a pos- 
itively charged D-particle and two baryons, The experimental search for the kind 
of the baryons may solve finally the problem of the isotopic structure of the D-particle. 
The two baryons being both nucleons it will indicate that the D-particles form a triplet; 
if one of the baryons is a nucleon, and the second a & (or A) hyperon, this will be an 
indication of the doublet structure of D-particles, and, finally, if one the baryons 
is a Z hyperon and the second a nucleon, this will be strong argument in favour 
of the singlet in the modified Rayski scheme. 

The search for supplementary support seems to be, at present, rather unneces- 
sary, since the experimental data are favouring the singlet assumption. 


V. Concluding remarks 


It is worthwhile to mention that the presented version of charge singlet yields 
the same practical consequences as the triplet in the scheme of Salam and Polking- 
horne [1]. Since the neutral D-particle in the former scheme has the same strangeness 
as well as other quantum numbers as the z?-meson, it should decay immediately 
into some z-mesons (the number of pions depends on the assumption on the scalar 
resp. ps-scalar character of the D-particle). The neutral D-meson would be thus 


unobservable in practice?). 


3) The distinction between our singlet and the Salam-Polkinghorne triplet would be, however, 
possible in principle, by means of a detailed analysis of the associated production of one nucleon and 
some pions in the pion-nucleon reactions. These reactions must proceed via a neutral D-meson (with 
strangeness 0) like the two examples below: 

a +p>n+D°>n+ka 

nt ل‎ n— p D? — p + ka 
It seems that the analysis of this possibility will be very difficult because of the correlations involved 
and the need to distinguish duch events from the production of pion showers. 
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ON THE THEORY OF SPIN WAVES IN FERRIMAGNETICS 
By Jerzy Kocısskı 


Ferromagnetic Laboratory, Physical Institute of the Polish Academy of Science, Poznan 
(Received March 14, 1960) 


We investigate a ferrimagnetic with spins pointing in the same direction, while 
at absolute zero temperature. As the model of this ferrimagnetic we take simple cubic 
and body centered cubic lattices with two kinds of spins s, S, attached to the lattice 
points of two magnetic sublattices in such a way, that a small spin s has as its nearest 
neighbours large spins S and vice versa. 

Spin waves are introduced in a manner analogous to that of Dyson for the ferro- 
magnetic case. As in Dyson’s theory, the Hamiltonian of the system is explicitly separated 
into a diagonal part representing the energy of free spin waves and a nondiagonal part 
determining the dynamical spin-wave interaction. A quadratic dispersion law with two 
frequency branches of spin wave propagation has been obtained. 


1. Introduction 


The subject of ferrimagnetism has been treated with quantum mechanical meth- 
ods by Vonsovskii and Seidov (1954), Van Kranendonk and Van Vleck (1958), 
and by T. Morita (1958). In the paper of Morita only an elaborate form of Hamiltonian 
for ferrimagnetic has been given. The first two authors have considered a crystal 
lattice consisting of three magnetic sublattices with spins s = 4, the spins in one 
sublattice pointing in the opposite direction than the spins in the remaining two 
sublattices. They have obtained a formula for frequency of spin-wave propagation, 
linear in dependence on the wave vector of the spin wave. This contradicts the results 
of Van Kranendonk and Van Vleck, and of semi-classical spin wave theory (Cofta 
1959). These authors have obtained quadratic law of spin wave propagation. 
The neutron diffraction experiments (Brockhouse 1957, Riste, Blinowski and Janik 
1959) speak in favour of the quadratic dispersion law. 

It is tempting to generalize the theory of spin waves in ferromagnetics in the 
form given by Dyson (1956), to the case of ferrimagnetism. We shall treat only the 
type of ferrimagnetic with two kinds of spins pointing in the same direction while 
at absolute zero temperature, and shall content ourselves for the present with separa- 
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ting the Hamiltonian into the diagonal and nondiagonal parts, and obtaining the 
dispersion law for the frequency of spin-wave propagation, which turns out to be 
a quadratic one. 


2. Definition of spin waves 
We shall discuss a finite, simple or body centered cubic crystal with periodic 
boundary conditions. Each atom is labelled by a lattice vector j and there is a fixed 
integer n such that the atom j + nk is identified with the atom j for any two lattice 


vectors j k. The total number of atoms assumed even is N = nš. 

We divide the crystal lattice into two magnetic sublattices in such a way, that 
an atom of one sublattice has as its nearest neighbours only atoms of the other one. 
We label the atoms of one sublattice with odd and of the other one with even integers. 


To each odd atom j there is attached a spin vector s and to each even atom j a spin 


vector S, both obeying the usual commutation rules which we write in the form 


[55 se] = Op s? [5 Se] = ó, S$ 
[s> sk] 0 O [ 2 x1 m 4 Sr (1) 
(st, sy] = 26% 5? [St, $7] = 26, 5° 
with 
RR, $t = $ + i$ (2) 


We take into account exchange interactions between nearest and next nearest neigh- 
bours. The interaction between next nearest neighbours is essential in obtaining 
dynamical interaction between spin waves analoguos to that of Dyson for the fer- 
romagnetic case. The Hamiltonian of the model has the form 


%=8+M (3) 
= -4 و[‎ [Xs Eget ES 5; o] = 
J 
Ebo S 544 P S Sa] (4) 
JA j4 
Mic ل‎ (5) 


In every bracket of (4) and (5) the first sum extends over odd and the second sum 
over even lattice vectors j* The sums of (4) extend also over all the vectors 6 which 


join an atom to its nearest neighbours, and over the vectors A which j join an atom 
to its next nearest neighbours. é represents an isotropic exhange coupling of magnitude 
ول‎ between nearest, and of magnitude J, between. next nearest neighbours. The sup- 
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position of positive J, and J, causes that the spins of both magnetic sublattices are 


pointing in the same direction while at low temperature. The magnetic energy M 


represents the effect of an uniform external magnetic field directed along the positive 
z-axis, of intensity 

s S 
H = — L= — L (6) 


Ms ms 


where m, ms stand for the magnetic moments of the corresponding spins. Magnetic 
dipole coupling between the spins is not included in this model. 
We have one reciprocal lattice for the two magnetic sublattices. It consists of N/2 


reciprocal lattice vectors defined as the set of vectors A such that for all lattice vectors 


j of the odd or even sublattices the quantities 


1 (7) 


n 

477 

are integers. We denote the reciprocal lattice vectors by Greek letters A, u, o, o, and 
A+n/2 u is identified with A. 

To every reciprocal lattice vector A there are attached two spin operators defined by 


REM LA 4 PI zw af NE 3; 28 
$4 = (i) >> exp (¿À J) Sj. S= B 3 exp (iA j) Sj (8) 
7 j 


E £ d d : 
with summation over odd j in the expression for s, and over even J in the expression 


for Se The commutation rules of the operators follow from (1): 


E NET Ware. 

DENE 8 Sin RE Sa | = (3) Saas 

^2 ^— N uS e a; N -4 A 

Ere $,] = — pl Satu (Si, Si |= — (3) Sith 

^ A i N S e دخ‎ N zt 2 
[Bts] = 2 (3) S [55-2 (3) 2m 9) 


The ground state of the whole system is the state |0 > defined by 


50-0 $F |0-—-—s|0» j odd 
$-|0)=0 S?|(0>=—S|0> j even (10) 
or equivalenty by 
N à 
3 020 Si )0( ع‎ — (3) s 00 |0> 
À NM 


for all A. 1 
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In analogy to Dyson's definiton in the case of ferromagnetism we define a spin- 


-wave state |1,» containing a single spin wave with wave vector À by: 
117» = (a st + C, 51) [05 (12) 


with numerical coefficients c}, C, depending on wave vector A. These states are orthog- 
onal to each other as well as to |», and may be normalised. In generalizing this 
expression to states with more than one spin wave we follow Dyson and define the 


spin wave state |a? containing a, spin waves with wave vector 7 by 


la) = ا‎ y NG; sf TC, S?)] |0> (13) 


For Xa,>1 the states (13) are neither normalised nor Vers to each other. 
Dyson's discussion of analogous states for the ferromagnetic case holds also in this 


case and will not be repeated here. 


3. Kinematical and dynamical spin-wave interactions 


The nonorthogonality of the states (13) produces an interaction between spin 
waves namely the kinematical interaction according to Dyson. 

Another spin-wave interaction arises form the fact that the Hamiltonian (3) 
is not diagonal in the states (13). This one has been called by Dyson the dynamical 
interaction. To obtain this interaction explicitly we shall calculate the effect of the 
Hamiltonian operating on (13). Using (8), the Hamiltonians (4), (5) may be written 


in the form | 
— Be TT > Gr aes) Saye a) 
= Xexp(ió-2, y =J exp (if -Ñ (15) 
3 risk Wis, 4 
M= L 0 (ê + SD | (16) 


To calculate Æ |a? we commute the operator Z through the operators (c; st + 


+C, $1) appearing in Eq. (13) until Z operates directly on the ground state DZ I 
This gives a sum of terms involving the commutators 


6, = 18, (ef + ارت‎ a + 6,8 — 
1 | 
de, in È Gn erd (Î êa — ê p+ m 


| cS Az 1 N -4 ا‎ 
+ (Ce Ya — Ce Ya—e) Asus Ten — „Ja (x) EN (Ya — ya-9) X 


4 
x [ey Sa — sf =) iC, (Si Stine Sa 
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plus a term E, |0 >, with E, given by 
H |0 > = E, |0 >, FE = — šJ, NSsy, — 4 Ja N (s* + S2) y, (18) 


Next we commute the Q, again through the operators (c, s + + C 0 until 0, 
operates directly on |0>. This gives a sum of terms involving the commutators 


= (077 (e sa C$?) = (19) 


Re 
1, (N =e € 
re ar X Dis Soa + Ta Sata Sha) — 


14 AE il DAR JO cross لمر‎ C, CES Eo Sau) 


with 
حارم‎ ONC) ON ale ee) (20) 
tiec SC es LIO ric Co) (21) 
T zm Ya En (fes umm Yrs SR Ve (22) 


plus a sum of terms (L + &,) |a» with e, given by 


Q, |05 = (L + e) (a 3 + €,57) |0) (23) 
Cy = + (s + S) [Ja Yo + Ja (yo — y] == + [(s — S)? [Jj yo — 
s J (Yo fir yo)? +4 Js Syi (24) 


calculated in the Appendix. 
The process of commutation comes to end after two steps because 


Rs (25) 


Now we shall discuss the term (19), which represents interactions between spin 
waves. For small wave vectors 0, o strictly in the limit for 0, o 5 tending to 0, this term 
is small, by virtue of the form of the coefficients I‘, I’, and Ti This is shown in the 
Appendix. Thus the interaction between long spin waves is MS as in the case of 
ferromagnetics. But the interaction term (19) differs from that for the ferromagnetic 
case in that it does not represent a simple scattering of two spin waves. The term‘ 
representing simple scattering two spin waves in the case of ferrimagnetics should 
have the form 


-1 ^ ^ ^ 
R,=- 5 5 ; 3 T, (6544 Sonat Cora Sia) (co Sr Coax Soa) (26) . 
e 2.427 3 - 
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with numerical coefficient I. The interaction term (19), contains the same operator 
components as (26), but with different numerical coefficients. The form of (19) may 
have some significance. For the present, we may replace the interaction term (19) 
by (26) with the coefficient I", so chosen that (26) is greater than (19). This is shown 
in the Appendix. š 
When the interaction. term (26) is used, instead of (19), the operation Z|a> 
may be written in analogy to the ferromagnetic case 


Flay = [Ey + )يه رج‎ + e] lay 22 041) (27) 


where the second sum extends over spin-wave states |b) which are obtained from 
la? by replacing one pair of spin waves (0, o) by a pair (0 — 71 o+ à). The Q, are 
numerical coefficients containing the T',. As in the case of ferromagnetics the Hamil- 
tonian (27) is explicitly separated into a diagonal part representing the energy of free 
spin waves and a nondiagonal part representing spin-wave interaction. This interac- 
tion is weak for spin waves with long wavelengths. The separation of Hamiltonian 
into diagonal and nondiagonal parts, does not depend on replacing of (19) by (26). 
The Hamiltonian remains explicitly separated also with use of (19), but it would 
then not be possible to write the interaction term in the form (27). - 

The formula (24) for energy e, of spin wave, represents the dispersion law, 
depending on the square of wave vector 0. Similar form of the dispersion law has 
been obtained for ferrimagnetics by means of semi-classical theory (Cofta 1959). 

The author wishes to express his thanks to Professor Szezepan Szczeniowski 
for helpful discussions concerning the subject of the paper. He is also indebted to 
Dr Henryk Cofta for comments on semi-classical spin-wave theory. 


Appendix 
A.Calculation of the energy e, of spin wave 
Operating with Q, — L(c, sf + C, S$) on |0) we obtain the expression 


([J (yg — Cole, * Ye) -- Ja s(yo— yle, 5? + [Jas(yo— cal C, y) + Ja S(yo — YC, 52)]0» 
EN 


The operation [Q,—L(c, St +C, S+)]|0> should yield the energy 7 Sf the spin wave 
(cosy + C, S7)|0 >; that means we should obtain an expression 


[Q, — L (e, $t + C, $10» — e, (e 3¢ + C, S*)| 0: 
This is possible only when the coefficients by Er $n and C, St in (Al) are equal: 


Ja S(yo — Cic, 3 Ye) H- Jas (yo — y) ia كول‎ (yo — Ca] Ge Š yo + C pre Yo), (A2) 
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This condition leads to a quadratic equation for c,/C, and yields 


مد .4588 --5—92( yl?‏ 9 حو SHS) [1 Go — 99 = vol V[n‏ دور 
Ja (A3)‏ 


with 
1 E Jai, 


After inserting (A3) into the coefficient of C, $t in (Al) we obtain formula (24) for 
the energy يع‎ of the spin wave. 


B. The magnitude of the interaction term (19) 


The coefficient I, may be written in the form 
T = Calo [Calco (va E Cole, i Va—o) z (Vira چ‎ Cole š الوك دور‎ (Bl) 


The ratios C,/c,, C,/c,, in the limit for p, o > 0, tend to 1 for the positive sign of the 
root in (A3), and tend to —s/S for the negative sign of the root. We cannot take the 
limit value —s/S since while inserted in (Al) it leads to finite positive energy for 
a spin wave with zero momentum. The limit value 1, leads to zero energy for zero 
momentum of spin wave. When in (Bl) the limit value 1 for C,/c,, C Je, is inserted. 
we obtain for 1: 


I = Ca Co (y; me OA VARIUS ADS Ens ) = Cç Co T (B2) 


where 17, is identical with the factor for the ferromagnetic case. It is small for small 
0, G irrespective of the value of A. The coefficients Cg, €, are determined by means 
of (A3) and the normalizing condition (1,|1;» = 1 of the spin-wave state |1,), 
and they are finite. Thus for small 0, o, the coefficient T} is small. The same holds 
for 15. 
The coefficient D is of the same form as I, and so is also small for small 0; 0. 
Thus the interaction term (19) is small for long wave-lengths of spin waves. 


C. The I" coefficient 


We may equate the coefficients of the corresponding operators in (19) and (26) 
to obtain the conditions to be fulfilled by the common factor I’,. We obtain 


T loa Sona == qnis. CAC (C1) 
TOR MGs E. COG, (C2). 
COE AA (C3) 


By bpp SS Tan, (C4) 
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These four equations yields different I", — s. Therefore (19) cannot be transformed 
to the form of (26). To replace (19) by (26) we may take as the common factor T, 
the sum of the four I, — coefficients determined by the four equations (Cl), (C2), 
(C3), (C4). This yields a term of form (26) larger than (19). 
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MOLECULAR THEORY OF LIGHT SCATTERING 
BY MULTI-COMPONENT SYSTEMS 


By S. KIELICH 
Institute of Physics, Polish Academy of Sciences, Poznan 
( Received March 31, 1960) 


A molecular theory of light scattering by multi-component systems of optically 
anisotropic molecules is proposed. General expressions, accounting for the degree of 
depolarization D of the light scattered, for Rayleigh’s ratio S, and for the turbidity h of 
the medium are derived. These expressions contain the molecular factors F;, and Finis 
of isotropic and anisotropic light scattering, respectively, which are discussed with 
special reference to multi-component systems of: (i) anisotropic molecules with constant 
polarizability, (ii) dipolar molecules with hyperpolarizability, and (iii) non-dipolar 
molecules with quadrupole moments. Formulas are derived relating anisotropic light 
scattering and the optico-optical birefringence of the medium, and others relating iso- 
tropic light scattering to the divergence between the molecular refraction of the substance 
in the condensed and gaseous states. Applications to one-component systems are given. 


l. Introduction 


Proceeding from ideas of Lord Rayleigh (1899) and Smoluchowski (1908), Ein- 
stein (1910) developed a phenomenological theory of light scattering by homogeneous 
and non-homogeneous liquids consisting of optically isotropic molecules. 

According to the Einstein-Smoluchowski theory, light scattering in solutions 
of liquids is due to local fluctuations of the density and fluctuations of the concentra- 
tion of the solution. The theory was verified experimentally by Raman and co-workers 
(1923). The phenomenological theory of light scattering by multi-component systems 
was further developed by Debye (1947), Kirkwood and Goldberg (1950), Stockmayer 
(1950), Fürth and Williams (1954), Tatsuo Ooi (1958) and Hill (1959). 

The theories of these authors dealt with isotropic light scattering only. In the 
general case, anisotropic light scattering, resulting from optical anisotropy of the 
molecules or from anisotropy of the molecular field existing within the medium, 
should also be accounted for. Thus, within a condensed medium consisting of aniso- 
tropic molecules, in addition to radial correlations, angular intermolecular correla- 
` tions intervene. This latter fact was taken into account by Benoit and Weill (1956, 


(575) 


574 S. Kielich 


1958) in computing the optical anisotropy of light scattered by dilute two-component 
systems of anisotropic molecules. 

In the present paper, a general statistical-molecular theory of light scattering 
by multi-component systems of optically anisotropie and polar molecules is proposed. 
The fundamental equation of the intensity of the light scattered is derived, containing 
the molecular factors F;, and F nis» which account for isotropic and anisotropic scat- 
tering, respectively. The equation is discussed for a number of special cases. 


2. The general theory 


Let us consider an isotropic medium as represented by a spherical specimen 
of macroscopic size. We shall consider the molar volume V therein, containing ولاليد‎ 
XN, ... molecules of components 1, 2, ... respectively, with N denoting Avogadro’s 
number, and x, — the mole fraction of the i-th component. The molar volume is 
given by 


TUM 2 x M;. (2.1) 


wherein summation extends over all the components of the system, M; denotes 
the molecular weight of the i-th component, and o — the density of the medium; 
the sum of all the x; is unity: 


» xc 


The theory of Rayleigh scattering by a system of the kind described, when 
in normal conditions, can be fully accounted for by the methods of classical electro- 
dynamics and statistical mechanics. Thus, we consider the system in the electric 
field E of an incident light beam of wavelength A. The intensity component of the 
light scattered by the volume V and passing through the analyzing Nicol prism at 
the point of observation is, quite generally, given by (see, Kielich, 1960): 


_ 1624 


B= = 
24 R 


«M, M; Na 000 (2.2) 
wherein R,is the distance of the point of observation from the centre of the scattering 
volume V, n — the unit vector perpendicular to the direction of observation (Rm = 0) 
and describing the plane of vibrations of the Nicol prism, whilst z and f are summation 
indexes assuming the values 1, 2, 3; the asterisk denotes the complex-conjugate 
quantity. 

M = M(t,E) is the dipole moment induced in volume V by the electric field E 
of the incident light beam, when the configuration of the system is v. In the classical 
treatment, the configurational variables + = t(r, œw) are continuous, and account for 
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the positions (r) and orientations (c) of all molecules in the system. For a multi- 
-component system, we have 


xN 
M (7, E) = 2; 2i me» ei* r(? (2.3) 


with m denoting the electric dipole moment induced by the electric field E of the 
incident light wave in the p-th molecule of species i, T? its radius vector, and, for 
Rayleigh scattering, 


> Ar . 0 
|s| = |k — k'| = Zinn (2.4) 


kand k’ being the wave vectors of the incident and scattered light, respectively, and 
V — the angle of scattering as subtended by the directions of propagation of the 
incident and scattered light waves. 

The symbol < >, in eq. (2.2) denotes the statistical average in the presence 
of the electric field E of the incident light wave. Let U(r, E) denote the total potential 
energy of the system in the configuration v and when the external field is E. By clas- 
sical statistical mechanics, the average statistical value of an arbitrary function of state 
Q(r, E) at thermodynamical equilibrium of the system is given by: 


U (t,E) 
f 96. E)e “kT dt 
TE er en a wert (2.5) 
سوط‎ doe | 


where the integrations with respect to ع‎ are over all configurations of all the molecules 
of the system, k is Boltzmann’s constant, and T — the absolute temperature. 

Restricting the problem to linear dependence of M(r, E) on the electric field 
strength E, we have 


9 M, s 
M, (x, E) = (Se), E,;. (2.6) 


accordingly, eq. (2.2) assumes the form 


16 n!I, /9M, [9 Mg" P 
bt RERO SEI ا‎ ee (2.7) 
where J, is the incident light intensity, and e — a unit vector having the direction of 
the electric field E. Here, the brackets ( ) with no index symbolize the average 


statistical value in the absence of the electric field (E — 0): 


ELOR 
tat YA 28) 


_ 0 (n0 
k 1 e RT dt 
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U(r, 0) is the total potential energy of the system when the configuration is + and 
E = 0. 

When the electric field E is absent, the probability for all directions of the unit 
vectors n and e is the same; thus, the product n,nge, e, of their components can 
be averaged isotropically, and eq. (2.7) yields 


16 1I, 
45 A* R2 
with Q,, denoting the angle between the vectors e and n. The quantities F;, and Fanig 
in eq. (2.9) are of the form 


T, = {5 cos? Qen Fis (s) + (cos? Ren + 3) Fania (s)} (2.9) 


xN xN 


8 m? [9 (4,3) * 5 x D 
AOE 3 2222 > ee ae) ى‎ D (2.10) 


i,j 


xN xjN (psi) (95) 
Oma Im, (pq) 
Fanis (5) = D) (a bay Ops — Sap o) 3 > 3E E s Sm Bi De VEM 


and are termed the molecular factors of isotropic and anisotropic light scattering, 
respectively ; r29 — r(9—1? is the vector connecting the centres of the p-th and g-th 
molecules of the i-th and j-th chemical species and 6,, — the substitution tensor 
(unity if x = f and zero when © = f). 

Eq. (2.9), together with (2.10) and (2.11), is the fundamental equation of the 
statistical-molecular theory of light scattering by a multi-component system. 

If natural light is used and the light scattered is observed with a Nicol analyzer, 
we have (see Kielich, 1960): 


16 z^ I, 
45 14 R? 


Io, = {5(1 — sin? g? sin? 9) Fis (s) + (7 — sin? g, sin? 9) Fanis (s)], 


(2.12) 


where رب‎ is the azimuth of the unit vector n, and  — the angle of scattering. 
If no Nicol prism is used, eq. (2.12) yields the following expression for the light 
intensity scattered at the angle 9: 


1(8) = {5 (E + cos? 8) Fis (s) + (13 + cos? 9) Fanis (s)). (2.13) 


45 = T 

The degree of depolarization is defined as the ratio of the smallest and largest 
possible values of the intensity of the light scattered. As Znin = 1(90°, 9) and I „= 
= I(0°, 9), by eq. (2.12) the general expression for the degree of depolarization 
of the light scattered at the angle 2 is obtained in the form 


5 cos? 9 Fis (s) + (6 + cos? 9) Fanis (s) 
ASS OE er 
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The light scattering constant, also known as Rayleigh’s ratio, is defined as the 


ratio of the intensity of light scattered by unit volume at angle 9, to the incident 
intensity: 


I (9) R2 
5 زق)‎ = VE (2.15) 


As the incident intensity existing within a medium of refractive index n is related 
to Í, by the formula 


3 2 
= f=) I (2.16) 
the following expression is obtained from (2.15) with respect to eq. (2.13): 


S (8) = eL. (5 (1 + cos? 9) Fis (s) + (13 + cos? 9) Fanis (s)}. (2.17) 


If the light scattered is observed perpendicularly to the direction of incidence, 
eqs. (2.14) and (2.17) yield: 


6 Fanis 


Dono ال‎ qo 
8 xš (n? + 2? | 
== — .x— H i is) 2 
S epo 0 Fis + 13 Fan). (2.19) 


Integrating S(9), as given by eq. (2.17). over the surface of the scattering sphere, 
the following general expression is obtained for the turbidity: 


2 
S K 1 (5 (1 + cos? 9) Fia (s) + (13 + cos? 8) Fania(s)}sin 849, (2.20) 
0 
or, in the absence of internal interferences: 
128 25 (n? + 2)? 
AN. DEO usse mno 2.21 
i 3i; y F] ) (2.21) 


In the subsequent sections, the factors Fj, and F,,;g will be discussed for several 
molecular models. 


3. Anisotropic molecules with constant polarizability 


In the first approximation, the polarizability of the molecules immersed in the 
medium may be assumed to remain unaffected by the surroundings, and thus to be 
that of the isolated molecules. In this case, if the molecules of each species present 
within the system under consideration are optically anisotropic, we have 
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with ag denoting the polarizability tensor of the p-th isolated molecule of species i. 
Hence, the molecular factors (2.10) and (2.11) assume the form 


XiN XjN 


s (ses > e» P at i) V cos S - voy. (3.2) 


Drake 


xiN x; N 
=1 > 2 3d (3 aor? ath = ap) o) coss: redy, (3.3) 


isj 1q=1 


FS 


Assuming all molecules of one species to possess the same polarizability, and 
applying the Kirkwood-Buff (1951) method, the molecular factor of isotropic light 
scattering may be written as follows: 


NA 
Fs (s) = A3 Ariy gy | 0; + x; x; n Tr, Ë (t) — » | d i} (3.4) 


bJ 


with a, denoting the mean polarizability of the isolated molecule of species i, and 
Zi, (z) — the distribution function of the correlation between the molecules of species 
i and j; 422 جل‎ = d r;; do; doy, where Q = f do. 

The molecular pio: ef anisotropic light scattering (3.3) will now be represented 
in a form better adapted to further discussion, namely: 


xiN xjN 
F nis (s) = 1 d 2 a) a S 3 mer ij) quei — bap ôa) coss: 2203 , (8.5) 
15J p=1 4-1 
with w%% denoting the cosine of the angle subtended by the axes of the molecular 
systems of reference.rigidly attached to the p-th and q-th molecules of species i and j, 
respectively. 
For molecules possessing the axial symmetry, this yields 


748 (s) = 9 N > Gi) Qj) Ost) Óp g 0; 5 + 
tJ 


22 = x, fe cos ?*0;; — 1) cos s ٠ z, ; au (z) — 2 a. (3.6) 


Here, ©,, is the angle between the axes of symmetry of the molecules of species 
i and j, respectively, and ó, — the optical anisotropy of the isolated molecule of 
species 7: 


(i) (i) G 
ô 5 aj Ee 8 a — aP (3 7) 
a 3%) E Hala ; 


% and ره‎ denoting the polarizabilities parallel and perpendicularly to the symmetry 
axis of the molecule, respectively. 
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Eqs. (3.4) and (3.6) can be rewritten as follows: 


Fi, (s) — 9 N Di xf xf ay as COR, (3.8) 
5j 
Be (sy oV. 2 xA x fj X, ریت‎ Sari) ریچ‎ ro (3.9) 


with the notation 
GE, * UH š > 
Cia = AM 4; 50; xP x f cosse Ti; | ريع‎ (7) — dt, (3.10) 


jy epe t 
CS à; af Besto- n ees ry |n (o 7 a8, 


(3.11) 


C62. and CEP c are quantities accounting for the radial and angular intermolecular 
correlations. 
If, in particular, a two-component system is considered, this yields 


NE. 
= Fig (s) = xoi 1 EET E Talen (7) — A aa + 


Nie 
+ 2 X1 Xo %(1) &2) fe S° Tio E (7) => d dt A 


A Nem 
+% 2) t + Xp E "Toe PË o= zl al, (3.12) 


eel ME ven xa Sait) fes + 3 1 (8 cos? 94, — 1) cos S: Ty, E (7) — A al = 


Ne 
€ Hy X &) %(2, Óa(1) Óa(2, f (8 cos? 0,5 — 1) cos S Tja E (z) = a dt +. 


a‏ را 
4g &(2) O22) h + > ii (3 cos? Oza — 1) cos S ` Tg, | Q= A a). (3.13)‏ + 


Hence, for a dilute two-component solution, and for A> Tijo the expression 
given by Benoit and Weill (1956) for the optical anisotropy 4? = F,nis/ is is derived, 
whereas, for a one-component system, the expression discussed by Benoit and Stock- 
mayer (1956) is obtained (cf. Anselm 1947; Buckingham and Stephen 1957; Kielich 
1958). 

In the case of an ideal gas, the radial and angular correlation factors (3.10) and 


(3.11) reduce to 


um at (3.14) 
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and eqs. (3.8) and (3.9) yield 
E ON > x; Oly = > x, FQ (3.15) 


gas 1s i is? 


SE ini = 9N > x; ; (59 ai) = = 22 x; pe (3.16) 


Thus, for an ideal gas, the molecular factors F;, and F,,;, result additively from 
the factors 


FP—9Nal and FY, = 9 Nat, 62.) (3.17) 


of the components of the multi-component system. This means that the various 
components of a gaseous system scatter light independently of one another. 


4, Molecules with permanent dipole moment and hyperpolarizability 


In their theory of the degree of depolarization, Buckingham and Stephen (1957) 
proved that in the case of a one-component dipolar liquid, the light scattered is largely 
dependent on the hyperporalizability of the molecules. Here, the method of the 
foregoing authors will be applied for computing the molecular factors F; and Finis 
of a multi-component system. 

In a condensed medium, the polarizability of a molecule is generally modified 
by the molecular field of its neighbours. In the case under consideration, the dipole 
moment induced by the total electric field (E+ F) in the p-th polar molecule of species 
i immersed in the medium is given by the expansion 


m»? = og ( y 4i TRD) +4 Pa ( TD as Ee) ( Ee 1: Fp) A (4.1) 


wherein Do is the er, tensor of the p-th molecule of species i, whereas 
F® is the molecular field at the centre of molecule p of species i due to the charge 
distribution of all the others. 


Substitution of the expansion (4.1) in eqs. (2.10) and (2.11) yields 


. SN N 
Fa) = DO es + a i + 
+ ag) Bab, FY + ...) cos s: POS, (4.2) 
XiN XjN j ; f 
F nis (s) = š$ P e 3 (a) d + a} BR FO + 
+ at BO. por Ue eget mn DE ag 2 cos S: reo’, (4.3) 


OF, 
The terms G, Sr have been omitted herein; these are discussed in Appendix A, 
g 


for the case of a one-component system. 


Molecular Theory of Light Scattering 581 


In the Onsager model (1936), the molecular field F® appearing in the foregoing 
expressions is the reaction field in the absence of an external field. The Onsager 
model was applied by Piekara (1950) to the theory of molecular orientation phenomena 
in two-component solutions of liquids consisting of anisotropic molecules. In the 
present case, the reaction field existing at the centre of a molecule of species i, of 
radius a;, is given as follows: 


G 
O 2E ال‎ ma O 
Fa = RZ = Sn = Mel, (4.4) 


wherein £ is the electric permittivity of the medium as a whole. In the case of dipolar 
molecules of anisotropic polarizability we have, at E = 0, 


m pus E) 1 ex DE (4.5) 


and eq. (4.4) yields 
; u 
Fo NT DD (4.6) 
1 — fia Ke 
wherein قر‎ is the component of the permanent electric dipole moment of an isolated 
molecule of species i and a(? — its polarizability referred to the principal axes. As- 
suming, according to Onsager, 


E. 


the quantity y® in eq, (4.6) can be written as follows: 


G) fita e 2(g — 1) (n? — 1) 


: ; aT 4.8 
a a NE BE nee) (n2 — T) a 8) 


where n, is the refractive index of the i-th component of the system, and A® = a. 
For molecules of isotropic polarizability, 4% = 1, whence eq. (4.8) reduces to 


£ : y  2(e— 1l) (2—1 ` 
Or DI AE LIE M ا‎ u 4.9 
Li 18: 3 (2e + nj) = 


By eq. (4.6), with eqs. (4.2) and (4.3), the following result is obtained for mole- 
cules whose polarizability has been referred to the principal axes: 


xiN N 


He > x44 2 coss: redy, (4.10) 


i,j æ, B=1 


xiN xj N 


3 
F . Du AG (M 2,0 cba ay (ps si) — ])coss: ae NS. (4.11) 
iri a p= 


EUN 
p=1 q-1 
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with the notation 


AG — a af SHC DBD, yD m — aa Bhoy Xy y ze = +... (412) 
1 


For light scattering by molecules possessing the axial symmetry, the foregoing 
expressions yield 


j. ] (ij) 
Fi, (s) = 9N > = x; xj (G) G) N + Xa £o Po bo sr YO) d Sr | (radiate 
i) (7) 
ij 


(4.13) 
y i) PG Ó ( 
Fanis (s) = 9N > xÉ % oa) Savi) Oat) 1 + XG) Hc) Êd Onc ES 
du) GC) Öali) 
yos ; 
+ xo APO PD 4, | NE (4.14) 
Qj) Oats) 
with xo = x9; 
(i) (i) 
E 1295. 
Osi) = Pi Ba and = By) = 3 (Bi) + 209 (4.15) 
3B) 


are the anisotropy and the mean hyperpolarizability of the dipole molecule. 

Thus, the effect of the hyperpolarizability and permanent electric dipoles of the 
molecules on light scattering by a multi-coniponent system has been accounted for. 
If, in particular, the scattering medium is a one-component system, the expressions 
derived by Buckingham and Stephen (1957)!) are obtained from eqs. (4.13) and (4.14) 
by neglecting the angular correlations and internal interferences: 


N 
Fa = 923 > n + 2% A RTE. (4.16) 


Fanis = 9a2 6 wh E ze (4.17) 


with R denoting the gas constant, مم‎ — the isothermal compressibility coefficient 
of the medium, and, for A; ~ 1, y given by 


IH 1) (n2 — 1) 
SR IDEE SUM A (4.18) 


1) Buckingham and Stephen (1957) computed the quantities A and B standing in connection with 
eqs. (4.16) and (4.17) at 
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3. Non-dipolar molecules with quadrupole moments and hyperpolarizability 


Let us now consider a multi-component system of non-dipolar molecules with 
hyperpolarizability. In the present case, the following expansion can be written: 


m» = icd E HORE RA LCS E: Fe») (EG) + pe?) 4 ar (E: is Ff?) 
(5.1) 
wherein y is the hyperpolarizability tensor of the p-th (non-dipolar) molecule 
of species i. 
By eq. (5.1), the molecular factors of isotropic and anisotropic light scattering 


assume the form 


xiN N : 
FG == = 2 P» 2 (at? aay +4 1 CS VEE NI (5.2) 


ij p-l 
+4 Pop +...) cos S: pv 
Fo) 41») or > {ag agg + + advan FO FU + (9.3) 
ij: p= q=1 
+4 et weed AUR URP + ...} ot99 quein — Jap Oya) COS S ٠ rio» . 
As previously, the molecular field F: can be replaced by the reaction field of 


Onsager's model. If, moreover, the molecule possesses an electric quadrupole moment, 
the reaction field at its centre is given by (see Appendix B) 


FO = f mf + h eros em 
wherein 
2(e—1) 1 6(e—1) 1 5 
a CN oe 


eo denotes the electric quadrupole moment of a molecule of species i, defined as 


follows: 
oO = 1 » ef? (3 r(9 ويم‎ 2). (5.6) 


where e? is the n-th ee charge of a molecule of species i, and TY? — its radius 
vector, summation extending over all the charges of the molecule. 

The reaction field of a permanent quadrupole induces a dipole moment in a non- 
-polar. molecule (the induced quadrupole moment will be neglected): 


| nf = a FP: | 57 
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FF oce 


thus, eq. (5.4) yields 


1 1 9% Ti 
FP = وير‎ = em SLM P (5.8) 
v Í 
wherein 
SEEN nal 
G) _ ads i) 1 (i) 
BT - (3225) (3) £e. 6.9) 


With the definition (see Böttcher, 1952) 


ai n 2n 


FO FO = = f I j FO FG 72 dr, sin ð; dd, dg, 
7 e e 


and with respect to eq. (5.8). the mean square value of the reaction field is obtained 


as follows: 


EG 1 nt 2| DD TN 
1 1 


For molecules having the axial symmetry, this yields 


Oy 


FO FO = 9,07 (5.11) 
as” 
wherein 
1 [25 E 1 : men 1 % (i)2 (2 


and © رم‎ is the quadrupole moment of the axially symmetric molecule, 
Oy = 2 4 (ef, — xin 


Zin being the distance along the molecular axis of the charge ef” and Xj, — its distance 
from the axis. 
With respect to eq. (5.11). eqs. (5.2) and (5.3) assume the form 


5 4 9. 
Fi (s) =9N ay x x Xei) Xi) h + > Nu) Oe s 
(i) 


ij 


5 yg0 à.) 
Toy 70 a. + of Cradials (5.13) 
G) 
I 7 sd 2. 
Fanis (s) = 9 N > xÉ xf Aci) Aj) dar) dal) 1 + + na TO RUE 9%) 5 
DT] 10) bari) 
1 yu ony OG) Gi) | 
Hec qe re | Cards 
4 Gs) c Ç E35) 
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which holds for axially symmetric molecules. Here, the quantities 
Hp eret 
3 yo» 
have been introduced, being the anisotropy of hyperpolarizability of the quadrupole 
molecule and its mean hyperpolarizability, respectively. 
Thus, the effect of the hyperpolarizability and of the permanent quadrupole 


moments of the molecules on light scattering has been computed. 


In the case of a one-component system, and neglecting the internal interferences, 
we have, by eqs. (5.13) and (5.14), 


Oni) = 


l oe à 
ME a: (5.15) 


N 9? 
Fy = 90th + 5 922 tos ET6r, (5.16) 
Lie o 7 760? | 
Fa c Satan T SPEO zu (5.17) 
where 
iY mq or = yt w 
dicum (5 pus 5) p T N (xi + 2 x9. (5.18) 
and | 
1 IV] : 
Cangular =1+ fe cos? O — 1) ES = 4 dt, (5.19) 


is the angular correlation factor (cf. Anselm 1947, Benoit and Stockmayer 1956, 
Kielich 1960). 


6. Relation between anisotropic light scattering and the effect of optico-optical 
birefringence in an isotropic medium 


From the discussion of sections 3—5, the molecular factors of isotropic and 
anisotropic light scattering are.seen generally to contain certain molecular constants 
accounting for the electro-optical properties of isolated molecules, such as the polariza- 
bility æ, hyperpolarizabilities 6 and y, and permanent dipole u or quadrupole © 
moment. The factor F; is, moreover, dependent on the radial intermolecular cor- 
relations, whereas F ni depends on the angular intermolecular correlations. For the 
greater part, the values of these molecular parameters are not known, and so the 
factors F,, and F,,,, are more conveniently expressed by quantities that are accessible 
to measurement. : 

Eq. (2.18) defining the degree of depolarization can be used for eliminating 
` the factor F, from eq. (2.19), which is that of the scattering constant; this yields 


20 8n (n? + 2)2 6 + 6 D 


' S= Fanis - (6.1) 


SPAS eS DEN 
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In considering light scattering by a multi-component system, it is convenient 
to introduce the molecular light scattering constant S,,: 


9V 
لا‎ a = ا‎ 6.2 
ins (n? a 2)? S$ ( ) 
by eq. (6.1), the latter may be rewritten as follows: 
8x4 6+6D 6D 
ek, == 9 9 24 iD.» Panis (6.3) 


A relation between anisotropic light scattering on the one hand, and the optical 
birefringence arising in an isotropic medium from the effect of a very intense light 
beam on the other, will now be derived. This latter effect was theoretically predicted 
by Buckingham (1956) and, as far as the present author is aware, has not been expe- 
rimentally detected as yet. Generalizing Buckingham's theory for multi-component 
systems, the molecular constant of optico-optical birefringence is obtained in the form 


P er mp 

Di G00 en 

m = qs (3 bay موة‎ — dan rs) p Qi ERE 34 
N xN 


(Psi) 9 (a) 
aL m2 (64) 


In deziving this expression it had been assumed that the frequency of the weak 
measuring light beam is identical with, or differs but little from that of the intense 
light beam giving rise to the optical birefringence of the medium. In this case, the 
molecular constant B) and the constant measured 


ال 
B? = 5 p (6.5)‏ 
are related as follows:‏ 
n?‏ 54 0 
B. = (n? + 2)4 VB, (6.6)‏ 


with 3? denoting the time mean square of E. 

In particular, neglecting the effect of the molecular field, the dipole moment 
induced in the p-th molecule of species i by the electric field E of the light wave 
can be expressed by the expansion (cf. Buckingham and Pople, 1955) 


ma = aap” Es + 1 BBY Ep E, + 7999 Es E, E, + ue (6:7) 
by eq. (6.4), this yields i 


o 2nN a 1 0 | 
By = ES b = وفع( أنه‎ + gm), w ويه‎ Oyo 6) Goa Mua). 


jj g=1 
(6.8) 
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Neglecting herein the term accounting for the hyperpolarizability y, the following 
expression is obtained for molecules presenting the axial symmetry : 


0 An N 3 
Bm = BET 2, ^9 «^ Óa(i) Óa( jy {Xi Ó; ; + 
17 
+ 1x,x; f (3 cos? 0;; — 1) g;; (z) dt); (6.9) 


for a one-component system, the latter expression becomes identical with the one 
derived by Buckingham (1956)?). 

In the case of scattering molecules whose linear dimensions are small as com- 
pared to the light wavelength A, the molecular factor of anisotropic light scattering 
as given by eq. (2.11) assumes the form 


XIN XN. (pA) N AD 
"E 1 d Omg dm, > 
Ins m 9 (3 Oy Ops Sr 0 Oye) Y 1 ) I 9E; ran . (6.10) 
tJ p=1 q— 


If, on the other hand, the first term relating to the hyperpolarizability of the 
molecules be omitted in the general expression (6.4) for B? , we have, by eq. (6.10). 


ee ENT Eon (6.11) 


Eq. (6.11) relates the molecular factor of anisotropic light scattering and the.molecular 
constant of optico-optical birefringence, and holds for an isotropic medium of arbit- 
rary density. 

With respect to eq. (6.11), the molecular constant of light scattering (6.3) is now 
given by 


1229 , .1+D 

Se EG KT B. (6.12) 

By eqs. (6.2) and (6.6), the latter relation can now be expressed in constants 
that are accessible to measurement: 

7223 n? kT (1 + D) po 

A4 (n? + 2? D 3 

The foregoing relations contain no molecular quantities explicitly, and hence 

can be conveniently subjected to direct experimental checking. If the values of S 

and D are known experimentally, eq. (6.13) yields that of the optico-optical birefrin- 

gence constant BP. On the other hand, for non-polar substances, B? can be expressed 


by Kerr's constant K: 


ge (6.13) 


(n? — 1) (n? + 2) K (6.14) 


ef BE TIERE 


2) The quantity L mu. by Buckingham (1956) differs from BP, by the factor 1/9, i. e. Be = OL, 
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7. Relation between light scattering and molecular refraction 
From eqs. (2.18) and (2.19). we have the following relation: 

s- 874 (n? + 2)? 6 + 6D 

ol Aye GOED 


wherein (6--6D)/(6—'7D) is the Cabannes factor (1929). In the foregoing relation, 
the factor F, is the only quantity to contain molecular parameters. It will now be 
proved that the molecular factor of isotropic light scattering F, can be expressed 
in terms of the molecular refraction of the medium 


Fis, (7.1) 


nè? — 1 


ll (7.2) 


For a multi-component condensed system, the molecular theory yields the 
following general expression for R,,: 


FN Die CE. (7.3) 


Hence, with the models discussed in sections 4 and 5, the following relations are 
derived for dipole molecules: 


Am pbu 
كار‎ E ai) + xo) HOPO 2 +. 3 
3 : |o Xa) Tm 1 (7.4) 


and, for quadrupole molecules: 


7 Oto 
Hiit — E > x; hots tis NG) (i) f 2 ze Z M (7.5) 


From these relations it is seen that the molecular refraction of a condensed 
system and the factor of isotropic light scattering given by eq. (4.13) or (5.13) are 
expressed in terms of the same molecular quantities. 


Now, restricting the problem, for simplicity, to that of a one-component system, 
we obtain by eqs. (7.4) and (7.5) 


Ry — gast up 
ZEN 3E Aisa 23 (7.6) 
RA LEE gaslt Ex 5 9? 
suis ndr LTE ra 
wherein 
47 
gas Po, = Y Na (7.8) 


denotes the molecular refraction of an ideal gas. 


Molecular Theory of Light Scattering 589 


The foregoing expression together with the factor F, as given by eq. (4.16) 
or (5.16) yield the required relationship: 


8 f Ree us „ RTBr 
Fis = 16 16 x° | Sr 2 | Angela) | it "NV a (7.9) 


which holds for either of the models under consideration. 


By eq. (7.9), the scattering constant of eq (7.1) can be represented by the follow- 
ing formula: 


N! 42 eS Bel, 


274 NV? 6570 s| za m RTT ` (7.10) 


lm‏ دقوع 

which contains no molecular quantities explicitly. 
Similarly, by eq. (7.9) and the relation between F nis and the anisotropic term 

in Kerr's constant, 


1215 n? (n? — 1) kTV 


Fanis = 2 z (e — ]) (e + 2) (n° +2)? 


we have, with respect to eq. (2.18): 
| و‎ [Bm — امالس‎ | _ _24 n? (n? — 1) V? Kanis 
gas Po, (e —1)(e + 2)(n? + 2)? eas R, 


The term (R, — gas) R, in the latter equation accounts for the effect of the 
molecular field on light scattering by a condensed medium. More exactly, terms of the 


(1.12) 


, ..., accounting for the anisotropy of the molecular field, have been 


neglected in the computation of Fj, and F nis in sections 4 and 5. In the first 
approximation, identical terms will appear in the expression for the molecular 
refraction R,. Divergences will arise only in higher order approximations, . but, 
` as these additional contributions are extremely small, they may well be neglected. 
Hence, the expressions of (7.10) and (7.12) still remain valid if the effect of anisotropy 
of the molecular field on light scattering is taken into account. In the present case, 
however, the following term due to anisotropy of the molecular field should be added 
to the expressions of eqs. (7.6) and (7.7), which account for the difference between 


the molecular refraction of a substance and that of the perfect gas: 
3 
NAG ee gas JU, 1 (a 25, 1 
i te ppp fous = Yar (7.13) 
| ER? 5 3a on reges a Ya 
wherein, for Onsager’s model, we have (see Appendix A) 


2 (n? — 1? (4, — 1) 
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For isotropically polarizable molecules A, = «,/x = l, and so the contribution 
resulting from anisotropy of the molecular field vanishes, i. e. y, = 0. 

Neglecting fluctuations of the molecular field and fluctuations of its anisotropy, 
we have 


A 
Rn. = = d (7.15) 


hence, eq. (7.10) yields the well-known Cabannes-King-Rocard formula (see Caban- 
nes, 1929): 
_ 7? (n* — 1)? 6 + 6D 


c HS STE 1 7.16 
> HN GE DAE C10 


whereas eq. (7.12) reduces to Gans' (1923) relation: 


6r D 242 n? Kanis 


6— 7D € (n IES 426402) (7.17) 


In considering a condensed medium consisting of spherical molecules. eq. (7.9) 


should be replaced by 


Re, Gradi (7.18)‏ | برهم وم 
V gas Pas‏ 


herein, we have approximately (cf. Buckingham and Stephen, 1957) 
Ry ig gasR,, 1 IF, & E 
ا‎ Rn ا کے‎ AR ee dT. 1.19 
ل‎ DE, Baer 


& (r) being the radial distribution function. 

Buckingham and Stephen (1957) proved that, in the case under consideration, 
the factor of anisotropic light scattering can also be approximated by the quantity 
of eq. (7.19): 


9 ES aap 
ps ——g? ده تحدم‎ 3 
5 N [ E ) (7.20) 


By substituting (7.18) and (7.20) in eqs. (2.18) and (2.19), the degree of depolarization?) 


°) Eq. (7.21) for the degree of depolarization D differs from Buckingham and Stephen's (1957) 
formula by the larger numerical coefficient (27 as against 7) of the second term in the denominator. 
The difference arises from the fact that Buckingham and Stephen neglected the effect of the 
molecular field on isotropic light scattering. If the term due to fluctuations of the molecular field is 
numerically negligible with respect to the one accounting for fluctuations of the density, i.e. if 
RE gas TU << RTT 


, eq. (7.21) reduces to 
gas! pi A 


f MIRE V ين‎ 
D=— m gas m Le uy m PR r“ g(r) dr. 
5 Rn RIG, CRIB 
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and Rayleigh's ratio are obtained as follows: 


6 | Rm -J CEDIA 
gas os 


(7.21) 
[10 == 2 IE Es 27 ps ` 


n‏ دقوع 


^s == aufm 
+ 33 pee ل‎ (0092) 


The foregoing expressions hold for condensed substances consisting of spherical 


molecules. For the ideal gas, D — 0, whereas eq. (7.22) reduces to Rayleigh's (1899) 
well-known formula 


_ m (n? + M jo RTór 
= 203 NV ruse 


CS a (n° — D* , _ 222 (n—1)? 


SN ener e cuu 


Finally, as Kerr's constant for a condensed medium of spherical molecules is given 
by 


_ an + 2)2 (e + 2)2 N ا‎ (Rm— gases 
Kee -Bls niy O كه‎ ED » (724) 


we have, by eq. (7.21), 


LIE d 27 qu P» Pr D |‏ اليه 


de 243 n? V l6z? NV (2 —9D)|' 


(7.25) 
wherein ,,,P,, = 4/3 Na, denotes the molecular electric polarisation of a gas con- 
sisting of spherical molecules of polarizability «,. 

If, for a given substance, Kerr’s constant K and the degree of depolarisation D 
are known experimentally, eq. (7.25) yields the mean hyperpolarizability y of a spher- 
ical molecule. 


8. Discussion and conclusions 


The formulas derived previously will now be applied to the aim of computing 
the hyperpolarizability of the carbon disulphide molecule. In the present case, we have 
the following numerical data: at t = 20°C, n = 1.636, e = 2.26, Bp = 92x 1071? 
CGS, D = 0.62 (Weill, 1958), and K = 12.1x 10-?? e. s. u. (Stuart and Volkmann, 
1933); by eq. (7.12), this yields 

Rm EX gas m 


— 0.097. 
gos Km 


By comparison with eqs. (7.) and (7.13) we have 
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and hence, for z = 8.77 x 10-?* cm? and the value of 7 = 0.007 as computed by 
3 
eq. (5.18) and š >) A, v, = 0.084, we obtain 
»-1 
yO? = 20.8x 10-55. 


On the other hand, the quadrupole moment of the CS, molecule can be computed 
numerically from the Debye formula (1920) for the constant a in the Van der Waals 
equation: 


wherein d is the diameter of the molecule. Since, for CS,, a = 11.3x10!? dyne 
cm! and d = 3.86 A, Debye's formula yields 


O = 28x 10-26 e. s. u. 


With the quadrupole moment as computed above, the mean hyperpolarizability of 
the CS, molecule is ; 


225591036 eXs dq. 


Measurements by Blaker, Badger and Gilman (1949) of the constant of light 
scattering in CS, yield a value of S = 151x 10-9 cm-1, (A = 4358 A, t = 25°C); 
thus, by eq. (7.10), we have 


Rm aes gas lim 
gast m 


On the other hand, with eqs. (7.7), (7.13) and the numerical values used above, this 
yields y = 36.5 x 107?6, a value almost 14 times larger than that previously computed 
from eq. (7.12). 

The quadrupole moment of the CS, molecule as evaluated above is certainly 
exaggerated. If a value of © = 6x 107? e. s. u. is assumed in place of © = 28x 10-6 
e. s. u., eq. (7.12) yields y = 56.5 x 10796 e, s. u. (instead of 2.6 x 10-39 e, s. u.), which 
is a value not far remote from that of y = 54.4 x 10-36 e, s. u. computed by Buckingham 
and Raab (1957) from Kerr's constant. 

For CCl, at t= 20°C, we have B = 105x 10-1? cgs, & = 10.5 x 10-?* cm}, 
n = 1.461, e = 2.24, D = 0.06 (Krishnan, 1925), K = 31.4 x 10-14 e. s. u. (Stuart 
and Volkmann, 1933); hence, eq. (7.25) yields a value of y = 4.7 10-36 e. s. u. 
for the mean hyperpolarizability of the CC1, molecule. Buckingham and Pople (1955) 
computed the mean hyperpolarizability of the CH, molecule from the molecular 
Kerr constant of the gas, obtaining a value of y = 2.6x 10 36 e, s. u. 

From the discussion of the foregoing examples it is seen that parallel theoretical 
and experimental investigation of light scattering in liquids can lead to the numerical 
evaluation of the hyperpolarizability of the molecules or of their quadrupole moments. 
The foregoing molecular theory of light scattering proves that investigation of light 


4 


— 0.265. 
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scattering in multi-component systems can provide equally interesting information 
on the electro-optical properties of molecules of various kinds and on the nature 
of the forces with which they interact. 

It results from the present paper that, in a multi-component system, in addition 
to light scattering on fluctuactions of the density and concentration and on those 
of the anisotropy and orientation of the molecules, an essential part is played by light 
scattering on fluctuations of the anisotropy of the molecular field and by scattering 
related to higher order effects such as hyperpolarizability of the dipoles or polariza- 
bility of the molecular quadrupoles. 
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Appendix A. The Effect of Molecular Field Anisotropy on Light Scattering 
in Liquids 
The terms omitted in the computations of sections 4 and 5 and relating to ani- 


sotropy of the molecular field will now be considered. The expansion of eq. (4.1) 
or (5.1), if higher order terms are neglected, yields 


Oma” _ a? 1 al goa FE 
OE, Zope 


whence eqs. (2.10) and (2.11) assume the form 


+. (A.1) 


xN xjN pO erae 
na = Yi CD aar s rss 
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QF 
The quantities — 2 Eg herein account for the anisotropy of the molecular field 
and are readily dud from Waage s (1936) model. Namely, we have 
m (A3 
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wherein 
o. 9 n° (nË + 2) — (n° + 2) [2 n° + 1) (nf + 2) — 2 (a° — 1) (ni — D AE] 
E (n? + 2) [(2 قم‎ + 1) (n? + 2) — 2 (n — 1) (n? — 1) 22] 


(A.4) 


is the parameter of the molecular field anisotropy. For a one-component system 
n, =n, and eq. (A.4) reduces to (7.14), whilst eq. (A.3) yields the expression (7.13) 


for molecular refraction. 
Restricting the problem, for simplicity, to the case of a one-component system, 
and neglecting short range angular correlatious, we have, by eq. (A.3), 


3 3 
2 E I N 
Fa = 9 a? mem 2 Bee + 3) —' > RTBr. 
g= 


@,B=1 


NE a — y Aa Àg + 2 (sve y hm) n 
«=1 «=1 


a, B—1 a, B—1 


3 
+3 2i y 


wherein y, is given by eq. (7.14). 


3 
J Mave val N, (A.5) 


a, B=1 


For molecules possessing the axial symmetry, the foregoing expressions reduce to 
aspa ?N 
Papa enc VIL SA TO r 


E Ay A š 
2 = 2 2 AV ier LL 


with Ay = As; Àp = 1, = 2, and Vi = Vs; VI = V. = ys. 


The effect of fluctuations of the molecular field anisotropy on isotropic light 
scattering (f;,) in benzene amounts to no more than 3 %; the respective figure for 
CS, attains, however, 18 %. Fluctuations of the molecular field anisotropy yield 
a far greater contribution to anisotropic light scattering (F nis), one amounting to e. g. 
28 % in benzene, and to as much as 87 % in CS,. 


Appendix B. On the Effect of the Molecular Polarizability Due to the Reaction 
Field Gradient on Light Scattering 


The author is indebted to Dr A. D. Buckingham for kindly directing his atten- - 
tion to the necessity, in the computations of Section 5, of accounting for the dipole 
moment induced in the molecule by the gradient of the reaction field (cf. eq. (5.4)). 
This involves the term containing the tensor B, 5s intervening in the effect of optical 
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birefringence induced in the medium by the gradient of a constant external electric 
field, as computed recently by Buckingham (1959). 

For simplicity, only a one-component medium will be considered. With Buck- 
ingham's (1959) general expression for the total dipole moment, m® the total diffe- 
rential polarizability of the molecule is obtained in the form 


= [s (p) a a. EP Te 3 i po Fe d 


1 OF 
E po " mlp) 8 i : 
er s Fys + Opy + —— Ben (B.1) 


The tensor B,,., is symmetrical in the pairs of indices «ff and y8, and accounts for the 
additional polarization of the molecule as induced by the field gradient F(? = OF US 
A discussion of the relations between the components of the tensors in eq. (B.1), 
for molecules having the spherical or axial symmetry, is to be found in the paper 
by Buckingham (1959). 

Substituting the expansion of eq. (B.1) in eqs. (2.10) and (2.11), and neglecting 
short range angular correlations, the following expressions for non-dipolar molecules 
are obtained: 


lA 2 = N 
TUS == 9 (88 si Xaa Y BByó Fy Es ES T Kaw Bas: yo Fys = 2 y RTBr, 
Pec =: i {3 Hep Xap 2 Ka & 9g s (3 KB Yaßyö ES X za V'8Byo) E. Fs gi (B.2) 
TOA HS 20 Bess) + Ab 


For molecules whose axis of symmetry is in the direction of the 3-axis, the fore- 


going expressions reduce to 


25 AT 
Fis = 3 æ {3 æ + (Yaaap 2 Yıras) Fa Fa d ej y RIBr. 
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Fais = (8j — i) (aj — & 4 (Vas — Yırap) Fe Br 
+ 5 وو8)‎ 8 — Bug) Fag + ++} N. (B.3) 
The former is identical with the one derived by Buckingham and Stephen (1957), 
whilst the latter one differs therefrom by an additional term containing the tensor 
B 


aß:yö 
"The field F, and field gradient F وي‎ are computed with respect to omen 


(1936) model: 
Ven RR (B.4) 
wherein (cf. Onsager, 1936) 
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is the reaction field, and (cf. Buckingham, 1959) 


6 (e — 1) Gap 
3e+2 a 


Rep = = h 0 (B.6) 
is the reaction field gradient inside a spherical cavity of radius a in a continuous 
medium of electric permittivity e due to a dipole and a quadrupole of, moments m, 
and @,, at the center. 

Neglecting the quadrupole induced, eq. (B.4) yields, for anisotropically polarized 
polar molecules, 


F, = Oup ra, (B.7) 


L "eni لل‎ 
Ve fede 
whence 


IF, h 


— > B.8 
Org 1 — fa Án Oas ( ) 


Fg 


Substituting the foregoing in eqs. (B.3) yields finally 


evt 9|1 157 E B y RTBr, (B.9) 
= 92 88 ig yog s o ! 
Pas = 9020, nee + SEE 31 : Š a de N, (B.10) 
where n is given by (5.18) and: 
1 2e+1\ IS 


Until the numerical values of the quantities 3 yû, = y, — y , and B = Bg; 
are known it is difficult to say which of the terms in eq. (B.10) plays the greater part. 
It seems reasonable that both are of the same order of magnitude. 
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Note added in proof. 


1. An account of the flucuation theory of isotropie scattering in binary fluid mixtures and of its application 
to critical opalescence is to be found in the following papers: Pearson, F.J. and Rushbrooke, G. S., 
Proc. Roy. Soc. [Edinburgh]. A 64, 305 (1957); and Pearson, F. J., Proc. Phys. Soc. 75, 633 (1960), 
which became available to the author thanks to the kindness of Dr. F.J. Pearson after the present 
had been handed in for printing. The foregoing papers contain i.a.a detailed discussion of the equation 
for I(s), as analogous to eq. (2.13), for the case when F;, is given by (3.12) and Fynis = 0. 

2. The computations leading from the general formulas (5.2) and (5.3) to eqs. (5.13) and (5.14) were 
carried out on certain simplifying assumptions not explicitly stated. These will be discussed in a sepa- 
rate supplementary note, wherein, moreover, the factors F;, and F,nis will be derived without these 


assumptions. 
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FADING OF TRACKS OF ALFA PARTICLES IN NUCLEAR 
EMULSIONS AND A METHOD FOR THE CORRECTION OF THE 
VISIBILITY OF ALFA PARTICLE TRACKS IN THE 
GAMMA-BACKGROUND 


Z. Lewanpowski, M. MAKOWSKA-RZESZUTKO, Z. WROBEL 


Institute of Nuclear Physics, Polish Academy of Sciences, Cracow, Laboratory for Nuclear Physics, 
Institute of Physics, Jagellonian ‚University, Cracow 


(Received August 29, 1960, translated paper received on Semptember 22, 1960) 


Photographic latent images of tracks of alfa particles and gamma radiation were 
removed in Agfa-K2 and HHK®H-K emulsions 100 y thick by using the vapour of 
hydrogen peroxide in different concentrations as oxidizing agent at 32°C and 28°C. 
The nuclear emulsions were saturated with a solution of thorium nitrate in order to 
obtain alfa stars. One part of each emulsion plate was irradiated additionally with 
gamma-rays from a 200 mg Ra source. Thorium decayed within a few days yielding 
stars. To fade the latent image of alfa particle tracks and gamma-radiation, the vapour 
of a 10% or 6% hydrogen peroxide solution at 32°C and 28°C for the HAK®H-K 
emulsion, and 1.5%, solution at 32°C for the Agfa-K2 emulsion, was used. 

Pieces of nuclear emulsion plates were placed in a closed glass vessel (exsiccator) 
above the H,O, solution for different time intervals. The temperature of the vessel was 
kept constant by means of the water jacket surrounding the vessel and connected to 
a thermostate. The plates, after being exposed and oxidized, were taken out and devel- 
oped simultaneously with a control plate which had not been oxidized. 

The dependence of the number of grains in the tracks (for alfa particles), or of 
the definite area (for gamma-rays), on the duration of vapour activity on the emulsion 
was investigated. A typical fading process was found: at the beginning of the oxida- 
tion process the fading rate of the alfa particle tracks is very small. It is only after 
a time interval exceeding half the time required for the total fading of tracks that 
the number of grains begins to decrease rapidly and the tracks very soon become 
invisible in the chemical background. A similar character of fading is obtained for 
gamma radiation; albeit, fading starts from the onset of oxidation and attains earlier 
than in the case of alfa particles the point from which it begins to decrease con- 
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Fading of Tracks of Alfa Particles 


602 Z. Lewandowski, M. Makowska-Rzeszutko, Z. Wröbel 


Because of the foregoing differences, the ratio of the squared number of grains 
per unit length of alfa particle tracks S to the density of grains of the gamma-ray back- 
ground T shows a strong increase in time, attaining its maximum near the point 
at which the grain density for gamma-ray drops to the grain density of the chemical 
background. The ratio S2/T can be used as measure of alfa particle visibility in the 
presence of gamma-ray background. For the improvement of this visibility, the present 
authors propose to oxidize an emulsion between exposition and development, e.g. in 
the vapour of a 6% hydrogen peroxide solution at 32°C for no longer than 55 minutes 
(during this time the number of grains of alfa particle tracks decreases but insigni- 
ficantly). 

For the oxidation conditions described above, there was no experimental exidence 
of the dependence of fading on the emulsion depth for the HHKDH-K emulsion. A gfa- 
K2 emulsion exhibited a weak dependence of fading on the depth (below 5% statistical 
error), but this only at the onset of the oxidation process (Fig. 1). 
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Fig. 5. The dependence of the ratio N/N, of grains in recoil proton tracks in HHKDH-K emulsion vs 
the duration of development, with N, denoting the number of grains in recoil proton tracks in 5 
-oxidized emulsion developed for 14 minutes. 

Graph „p“ is for a part of the emulsion near the surface. 

Graph ,c* is for a central part of the emulsion. 
Graph „s“ is for a part of the emulsion near the glass. 


Typical experimental fading graphs are shown in Figs. 1, 2, 4 and 4. 

It would be of interest to test how the sensitivity of the emulsion was modified. 
The Agfa-K2 emulsion, after 7 hours oxidation with 3% H,O,, was able to register 
protons. In order to investigate the decrease of sensitivity m 1 DI]-K emulsion, 
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these were oxidized in the vapour of a 10%, solution of hydrogen peroxide at 32°C 
during 80 minutes, and subsequently exposed to Po+Be neutrons. The number 
of grains in recoil proton tracks was measured in 60 بر‎ length at 60 u residual range. 
Next, this number was compared with the analogous number in non-oxidized emulsion 
developed for 14 minutes. The obtained weak dependence of the sensitivity on the 
depth of the emulsion did not exceed statistical error. A decrease of sensitivity of 
about 20% due to the oxidation was found. This decrease can be reduced by developing 
for longer time intervals (Fig. 5). 

The authors succeeded in obtaining a method of fading alfa particle tracks, for 
HHKQoH-K and Agfa-K2 emulsions, in a time of less than two hours, without decrease 
of sensitivity. The method proposed is much faster than those known hitherto and 
can be very useful for work with old emulsions with high radioactive contamination. 

Moreover, the authors propose a method for increasing the visibility of heavy 
particle tracks in the presence of the gamma-ray background. 

A more detailed account will be published shortly. 
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